Mon. Not. R. Astron. Soc. 000, 1-30 (2009) 



Printed 4 April 2009 



(MN WT&i style file v2.2) 



On the variation of black hole accretion disc radii as a 
function of state and accretion rate 

O C. Cabanac 1 *, R. P. Fender 1 , R. J. H. Dunn 1 ' 2 ! and E. G. Kording 1 ' 3 

f"*"^ , School of Physics and Astronomy, University of Southampton, Southampton S017 IB J, UK 
■ 2 Excellence Cluster Universe, Technische Universitdt Miinchen, 85748 Garching, Germany 

3 AIM - Unite Mixte de Recherche CEA - CNRS - Universite Paris VII - UMR 7158, CEA-Saclay, Service d'Astrophysique, 
F-91191 Gif-sur-Yvette Cedex, France 

' Accepted 2009 April 02. Received 2009 April 02; in original form 2008 November 27 



ABSTRACT 

In response to major changes in the mass accretion rate within the inner accretion 
(-h ' flow, Black hole binary transients undergo dramatic evolution in their X-ray timing 

O !■ and spectral behaviour during outbursts. In recent years a paradigm has arisen in 

which 'soft' X-ray states are associated with an inner disc radius at, or very close 
to, the innermost stable circular orbit (ISCO) around the black hole, while in 'hard' 
X-ray states the inner edge of the disc is further from the black hole. Models of 
advective flows suggest that as the X-ray luminosity drops in hard states, the inner 
disc progressively recedes, from a few gravitational radii (R g ) at the ISCO, to hundreds 
of R g . Recent observations which show broad iron line detections and estimates of the 
disc component strength suggest that a non-recessed disc could still be present in 
bright hard states. In this study we present a comprehensive analysis of the spectral 
components associated with the inner disc, utilising bright states data from X-ray 
missions with sensitive low-energy responses(e.g. Swift, SAX), including reanalyses 
, of previously published results. A key component of the study is to fully estimate 

systematic uncertainties associated with such spectral fits. In particular we investigate 
\J | in detail the effect on the measured disc flux and radius of having a hydrogen column 

density that is fixed or free to vary. We conclude that at X-ray luminosities above 
~ 0.01 of the Eddington limit, systematic uncertainties only allow us to constrain 
the disc to be < 10i? g from spectral fits. There is, however, clear evidence that at 
X-ray luminosities between 10 -2 -10~ 3 of the Eddington rate, the disc does begin to 
k>( \ recede. We include measurements of disc radii in two quiescent black hole binaries 

. at bolomctric luminosities of < 10~ 7 Eddington, and present the inferred evolution 

C$ ' of disc luminosity, temperature, inner radius and accretion rate/efficiency across the 

entire range of bolomctric luminosities 10 _8 -1 Eddington. We compare our results with 
theoretical models, and note that the implied rate of disc recession with luminosity is 
consistent with recent empirical results on the X-ray timing behaviour of black holes 
of all masses. 
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1 INTRODUCTION 

The large changes in the X-ray spectra, radio fluxes and 
timing properties observed among the states of the BHBs 
are commonly explained with models involving a variation 
of the accretion efficiency and rate. This variation of the 
efficiency can be caused by the inward or outward motion of 
the inner radius of the optically thick accretion disc i?j n (see 
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e.g.. lEsinl Il997h . Recently, a number of papers have been 
published contesting this radius evolution. 

Since the first X-r ay spectral stu dies of Cyg X- l 
(|Tananbaum et al.l [1973 ) and A0620-00 (|Coe et all Il976t ). 
BHBs in outburst seem to transit between what are now 
commonly identified as "canonical states" . Those variations 
were subsequently associated with changes in radio emis- 
sion and the timing characteristics of the X-ray emission. 
A dichotomy could clearly be identified, resulting in the 
the characterisation of two states. In one state the spec- 
trum is dominated by a thermal optically thick component, 
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is hereafter called the soft or " Thermal Dominant" (TD) 
(|McClintock fc Remillardll2003rJ ) state. The other state, the 
so-called hard state, has a powerlaw shape X-ray spec- 
trum extending up to a few hundreds of keV (see e.g., 
iMcConnell et ahll2002l for the differences observed in Cyg 
X-l spectra). 

Those objects exhibit ejections that are observed in 
radio and infrared and that occur in particular states 
l|Fender fc Bellorn] 12004 ). Although steady and compact 
jets are typical of the hard state, the ejection process 
appea rs to be quenched whenever it enters in t h e sof t 
state (|Fender et all 119991 ). iFender. Belloni fc Gallcl (|2004T l 
give simple physical interpretations of the observed corre- 
lations, resulting from which promising elaborate models 
based on accretion-ejection s olutions in magneti sed discs 
have been built-up (see e.g., iFerreira et al.l I2OO6J and ref- 
erences therein for an auto-similar analysis of the problem, 
but also iMachida. Nakamura fc Matsumotoll2006h . 

For the vast majority of models invoked, the transition 
between the two main canonical states has been interpreted 
as follows. In the soft state, usually occurring at a high lu- 
minosity (i.e. 0.01 < -Lboi/-I'Edd < 1), the geometrically thin 
and optically thick accretion disc is believed to reach the In- 
nermost Stable Circular OrbitQ (ISCO). However, in order to 
reproduce the decreasing disc efficiency observed in the hard 
states, the inner part of the disc is truncated and replaced 
by a radiatively inefficient and optically thin flow. The ac- 
creted mass is either adv ected towa r ds the central black-h ole 
(ADAF and equivalent: lEsinl 1 19971 ; iNaravan et al.lll997l ) or 
part of it is ejected in outfl ows, for example, an ADIOS 
jBlandford fc Begelmanll2004 ) or a Jet Emitting Disc (JED, 
IFerreira et alj [20061 ) are examples of structures that could 
explain the observed behaviour. 

XMM-Newton observations of GX 339-4 |Miller et all 
l2006al . hereafter M06) cast doubt on this accepted interpre- 
tation by showing that a broad iron line together with a 
dim, hot thermal component was present in its spectra dur- 
ing the hard state. This effect seems to be observed in a few 
other sources such as Cyg nus X-l, and SWIFT J1753. 5-0127 
(M06. lMiller et al.|[2^06bh . However, all of these conclusions 
are based on single observation without studying the overall 
evolution during an outburst of the source. 

A robust way of evaluating the disc inner geometry 
value is to choose the same approach as lGierliriski fc Done. I 
|2004l ). They followed the disc parameter value during the 
outburst rise and decline. In the simplistic hypothesis of an 
optically thick disc extending down to the last stable orbit 
of the black-hole, one would expect that the relationship be- 
tween the disc luminosity Ldisc and the inner temperature 
Tin should be monotonic; Ldisc = K^ sc T-^ n with iodise con- 
stant, whatever the accretion rate. When the value of Ri n is 
much lower th an -R out , then i odise ~ 4mB% n <TB- 

Recently, iRvkoff et all (|2007l ) (hereafter R07) per- 
formed a similar analysis on the black hole candidate XTE 
J1817-330. They reported that no noticeable changes could 
be observed in the inner disc radii values, even when reach- 
ing the hard state. O n the contrary, taking int o acco unt the 
effect of irradiation, iGierliriski. Done fc Page! (|2008l ) argue 



1 The ISCO isequal to 6 R g for a Schwarzschild BH, or 2 R g for 
a fully spinning Kerr BH. 



Table 1. Orbital parameters used for the whole study. 



Source 



Mass (Mq) Distance (kpc) 



i n 



J1118+480 


(8.53±0.6) a 


(1.72±0.10) a 


(68 ± 2) a 


GX 339-4 




(8§;f) c 


(30|°) d 


J1655-40 


(6.3±0.5) e 


(3.2?;|)/ 




J1817-330 




(6.3^° 6 )' 1 


(60) h 


J1753.5-0127 




(2.9^,)* 




A0620-00 


(11 ± 1.9)-» 


(1.16 ±0.11)-» 


(40.75 ± 3)-» 



Gelino et ail (|2006h . 



Hvnes et all ||2003|) for the mass function. 
Zdziarski et al.l ||2004| ). 



b 

c 

d Compatible withlCowlev et al l l|2002h . 

e iGreene. Bailvn fc Orosj ll200ll ) 

/ iHiellming fc Ruper] l ll995h for the value but we used lFoellmi et all 

as a lower limit . 
9 Ivan der Hooft et al.1 lll998l ). 

^ S07 for the mass and distance, however see the text for the lower 
limit in distance. Arbitrary inclination. 

1 Mostly arbitrary but average values expected for BHB. 
Mill er et all ||20 06b ) for the low limit in distance. 
iGelino et aD l|200lh . 
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that a significant increase in R[ n could be measured. One 
of the goals of this paper is to disentangle this apparent 
discrepancy. 



Source sample and data selection criteria. 

Considering the large amount of data available in archives, 
it was necessary to select the data used in this study. The 
data were selected according to the following criteria: first, 
we focused on X-ray transients, as they allow to follow a 
source on broader luminosity range. As we focus on the low- 
est luminosity states and the study of the thermal compo- 
nent, data coming from instruments with good sensitivity 
at low energy (typically under 1 keV) were preferred, ruling 
out e.g. RXTE archives. For these reasons, the Swift/XRT 
data seemed to be ideal as it tends to follow the outburst 
of a source by taking a number of different snapshots over 
time. A major part of this analysis will therefore be based 
on data from this instrument. Lastly, confirmed black hole 
binaries were favoured, as we wanted to make comparisons 
among different sources using e.g. Eddington luminosities 
and hence masses and distances. However, as the outbursts 
of XTE J1817-330 and Swift J1753.5-0127 were well sam- 
pled, we also added them into our sample even though they 
have not been confirmed as black hole binaries. According 
to these criteria, our sample reduces to 6 sources: XTE 
J1118+480, GX 339-4, GRO J1655-40, XTE J1817- 
330, Swift J1753. 5-0127 and A 0620-00 (see Tab.[Qfor 
the adopted masses, angles and distances). 



2 DATA REDUCTION AND REANALYSIS 

The purpose of our paper is to extend previous studies of 
the disc parameters to lower disc luminosities, as well as 
to systematically and uniformly test spectral models. There 
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are two approaches that can be chosen to overcome the diffi- 
culties of determining the parameter values of the disc. The 
first is to try to fit the spectra with the best models avail- 
able. This means using sophisticated comptonisation codes 
or tables to fit the hard component, and the possible repro- 
cessing of hard X-rays on the optically thick disc. The main 
advantage is to usually obtain better fits to the spectra. This 
comes at a price, of having a larger number of free parame- 
ters. The other main drawback when adopting this approach 
is that comparison among sources or different observations 
is difficult as the models chosen are not necessarily the same. 

The other method consists of using the simplest model, 
involving the least number of parameters. This allows the 
same simple model to be fitted to the spectra resulting from 
every source and observation. The price which is paid is 
a lower goodness of fit, but it allows comparisons among 
sources to be made. 

We choose to adopt this latter approach in our study. 
The spectra were fitted with a powerlaw at high energy and 
any possible soft component by a multicolour disc. We also 
added the photoelectric absorption of the neutral hydrogen 
in the line of sight. Therefore the number of free parameters 
in the model does not exceed five: the absorption Ah, the 
multicolour disc temperature T ln and normalisation iodise, 
and the powerlaw photon index F and normalisation -Kpl- 
However, in order to study the effect of changing the high en- 
ergy emission model on the conclusions drawn, we also used 
a thermal comptonisation model (compTT in xspec) for some 
of the observations. As the value of the high-energy cutoff 
was usually unavailable, we therefore fixed the value of the 
temperature of the thermal electron to 50 keV unless stated 
otherwise. The number of free parameters are therefore still 
reduced to two for this component: the optical depth and 
the normalisation. 



2.1 Precision of the cross-analysis 

In this paper we pay particular attention to comparisons 
among several sources (Section [SJ), which involves differ- 
ent orbital parameters and distances, disc radius or lu- 
minosities have all been rescaled in terms of the gravita- 
tional radius R g = GM/c 2 = 1.48M/M km and Edding- 
ton luminosity L Edd = 1 .48 X 1O 38 (M/M ) erg s" 1 (e.g., 
iGierlifiski fc Done. Iliool ). Depending on the disc inclina- 
tion and for projection reasons we use: 



Ldis 



F disc 2nd 2 

X 



cos(i) 
4 x 10 5 x Fdisc dkpc 



Mm cos(i) ' 



(1) 
(->) 



where dk pc is the distance expressed in kpc, -Fdisc the disc 
luminosity in erg cm~ 2 s -1 , Mm q the black hole mass ex- 
pressed in solar masses, and i the inclination angle (i = is 
a "face-on" disc). 

Assuming that the hard component emits isotropically, 
we have approximately: 



ihard „ , „5 ^hard ^kpc 

— = 8 x 10 x 



(3) 



model (diskbb in xspec, see e.g.. lMitsuda et al]|l984l '). Sev- 
eral drawbacks have already been pointed out concerning 
this model as it neglects processes such as increased scatter- 
ing at high temperatures, or the differential emission that 
can occur between different altitudes in the disc (the central 
part of the disc may be warmer). Such drawbacks mainly 
affect the effective temperature compared to the observed 
one. However, the temperature shift is only 4% for a 0.5 
keV disc and this effect decreas es with decreasing tempera- 
ture (|Gierliriski fc Done. Iliooih . The fact that we aimed to 
probe the lowest luminous states also motivated our choice 
to use a simple diskbb to model the thermal component. 

As Rin is obtained via the normalisation (K d i BC = 
(Rin, km/dio kpc) 2 cos(i) in xspec) of the multicolour disc 
model, it gives: 



Rh 



= 0.677 



R g V 10k P c 7 V M ©7 \cos(i) 



M 



-1 / K diB , 



(4) 



For the evaluation of the disc parameters, and for the rea- 
sons explained above, we use the simple multicolour disc 



Unless mentioned explicitly in the text, all errors in this pa- 
per will also be expressed in terms of 90% (1.64 a) confidence 
range (including the plots). 



2.2 XTE J1118+480: multiwavelength campaigns 

2.2.1 In outburst 

For the study of the BHB XTE J1118+480, we mainly used 
the parameters given in the literature. The number of ob- 
servations where a thermal component is required in the 
spectra are quite scarce, as there is just the one coming 
from the 2000 outburst and one other, during quiescence in 
2002. Among the articles dealing with this outburst where a 
thermal component is observed, two were based on a multi- 
wave leng th study in c ludin g UV data l|McClintock et all 
l200ll and lChatv et all (|2003l ). hereafter (MC01) and (Ch03)) 
where as the other focused on results given by Beppo-SAX 
data (jFrontera et al]|2003h . 

We note a particular issue regarding the parameter of 
the disc obtained by MC01. There seems to be an inconsis- 
tency between the low value of the radius given in the text 
when using a simple multicolour disc to fit the HST and 
EUVE data, and the spectrum given in their Fig. 3. They 
obtain an internal radius of 34 Rs (IRs = 2R g ), and an 
internal temperature of 24 eV. Such a model, (assuming as 
(MC01) do, M = 6M , d = 1.8 kpc and i = 80°) should 
peak at a value of vF v ~ 10~ n ' 8 erg cm -2 s . From their 
Spectral Energy Distribution (SED) it is clear that this is 
about two orders of magnitude lower than the peak value 
shown in the third panel (i.e. vF v ~ 10~ 8,6 erg cm~ 2 s _1 ). 
In order to obtain a similar value with their orbital param- 
eters and temperature, we have to set Adisc ~ 7 x 10 8 . This 
then gives us R[ n ~ 300i?s when we take into account the 
orbital values displayed in our Tab. [T] 

This is moreover consistent with the value of R[ n ~ 
352i?s obtained by (Ch03) based on the same data set 
(HST, EUVE, CXC, RXTE, UKIRT and Ryle telescope), 
where SAX and the VLA data were added, and modelled by 
the same simple multicolour disc + powerlaw(s). The small 
residual differences could come from the value of the ab- 
sorption adopted by each authors (Ah = 1.1 X 10 20 cm -2 
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for (Ch03) and Ah = 1.3 x 10 20 cm" 2 for (MC01)), and the 
secondary contribution in optical as well. 

However, when using SAX data in their multi- 
wavelength analysis, (Ch03) did not take into account the 
lowest energy bins (0.13-26 keV) as the flux value conflicts 
with the EUVE ones. They argued that this discrepancy 
could come from calibration issues in SAX responses at low 
energy. However, as the amount of available observations 
for XTE J 11 18+480 outburst is quite low, we have tried to 
include these data in our study. This is also supported by 
the fact t hat in studies based on th e full SAX broad band 
spectrum (|Frontera et al. 2001, 20Q3|), a thermal component 
was also detected though peaking at a higher energy (~50 
keV) than tho s e give n in (MC01) and (Ch03). However in 
IFrontera" et all (|2003h . the low energy part of the spectrum 
was fitted with a single black-body spectrum and thus we 
tried to probe the effect of reanalysing the data with a mul- 
ticolour disc. 

Therefore, we used and reanalysed the processed 
LECS (low-energy concentrator spectrometer, 0.1-4 keV, 
iParmar et all 119971 ). MECS (medium -energy c oncentrator 
spectrometer, 1.8-10.5 keV iBoella et afl 1 19971) and PDS 
(Pho swich Detector System, 10-200 keV, Frontera et al.l 
1 19971 ) spectra available on the archive for both "ToOl" (Ob- 
sld 21173001) and "To03" (Obsld 2 11730012) observations 
mentioned in IFrontera etafl (|2003l ). The usual response 
files were used, we allowed a free normalisation among the 
three instruments and on e per cent syst e matic s were also 
added. As pointed out by IFrontera et all (|2003l ). the addi- 
tion of a thermal component to the absorbed cutoff power- 
law model is necessary. When we add a diskbb (the ab- 
sorption being free to vary), the reduced x drops from 
X jv = 286.05/191 = 1.50 to 213.64/189=1.130 (F-test 
probability = 1.05 x 10~ 12 ) in ToOl (from ^ /v = 2.98 
to 1.25 in To03) . We noted also that taking a single black 
body for the thermal component gives similar goodness of 
fit as in that case X 2 /v = 212.60/189 = 1.126 in ToOl for 
example. As stated before, we thus kept the results coming 
from a multicolour disc. The results of the fits to the disc 
parameters are given in Tab. [C]of the appendix. 

We note that the results we obtain for the thermal com- 
ponent are quite different between the SAX data and the 
EUVE data, as was pointed out in (Ch03) . The value of the 
disc flux is about 10 times lower in SAX data, the temper- 
ature is about three times higher than in the EUVE. How- 
ever the major difference comes from the normalisation of 
the multicolour disc as it is 10 3 lower. This means the disc 
is closer to the central black hole by a factor 30. As there is 
no obvious reasons to favour one or the other set of data, we 
decided to keep both of them, and infer that the real state 
of the disc should be between both. In the following figures, 
results coming from both data sets are displayed explicitly. 



2.2.2 In quiescence 

For the analysis in quiescence, we mainly based our anal - 
ysis on the results coming from IMcClintock et al.l (|2003al ) 
(MC03a). Compared to those authors, we did some slight 
corrections to estimate disc radii and luminosities for XTE 
J1118+480, given th e accuracy of the or bital parameters ob- 
tained since then bv lGelino et al.l (|2006l ). For example, with 
the value of the internal radius that (MC03a) obtain (around 



3100 R g but for M = 7M Q , d = 1.8 kpc and i = 80°), we 
derived the corresponding diskbb normalisation (.Kdisc). 

We then added the temperature estimate to evaluate 
the net disc flux contribution (Fdi sc oc itdiscT 4 ). Subse- 
quently we corrected the luminosities and the internal ra- 
dius from the newest orbital parameter in order to obtain 
R in ~ 1640 R g , Lhard ~ 1.3 x 10" 8 L E dd and L hol ~ 
7.7 x 10~ 8 I/Edd- Note that, except for the disc tempera- 
ture, no statistical errors coming from the fits were available 
for the luminosities nor the radius inferred and hence those 
latter are not displayed in the figures of the study. 



2.3 A0620-00: HST/STIS and Chandra 

In order to estimate the disc properties of A0620-00, we 
used the v alues published by (MC 03a), IMcClintock et all 
|2000l ) and IMcClintock et all i| 19951 ) as following. We have 
already pointed out that the lack of optical-UV observa- 
tion of the source (compared to XTE J1118+480) is still 
a major problem as we are not sure that a multicolour 
disc is the best model to fit the data. We noted that 
MC03a, IMcClintock et all (|2000l ) interpreted those data in 
the framework of an ADAF, whereas in our study we try to 
explain the possible UV extra component in the context of 
a multicolour disc. 

Another limitation concerning this source comes from 
the fact that the X-ray and optical observations were not 
simultaneous. Therefore the conclusions resulting from this 
object must be taken with caution. 

However, concerning the disc properties, and due to 
the simplicity of the model, we only had to find a rough 
estimate of the inner temperature and the normalisation 
of the diskbb. For the temperature, we note that the 
HST/ STIS disc spectrum obtained by IMcClintock et al] 
|2000l ) in 1998 was very sim ilar to the one obtained six years 
befo re (comparing Fig. 2 of | McClintock et al.ll2000l and Fig. 
5 of IMcClintock et al.l H~995l ) . Hence, we took the 9000 K 
(7.7 x 10" 4 keV) black body spectrum |McClintock et al.1 
Il995l ) as a reference for the temperature value. A diskbb 
spectrum with such a temperature value peaks nearby 3500 
A, as it was also observed in the 1998 spectrum. 

This temperature value fixed, we then obtained the 
normalisation of the diskbb model via the flux value at 
4 x 10~ 3 keV, which is around vF v ~ 10~ 12 ' 4 erg cm~ 2 s _1 . 
ft thus gave us -ft'disc ~ 2.8 x 10 11 . Hence, the corresponding 
radius would be (following mass, distance and angle value 



of Tab. [[]): R in = 4.5 x 10 3 R g . The bolometric flux of the 
disc component is i*di sc = 2.3 x 10 -12 erg cm -2 s _1 . 

For the hard X-ray component, we took the powerlaw 
model given by (MC03a) where the photon index value is 
given as V = 2.26. We inferred the normalisation value 
of the powerlaw considering that the flux value at 1 keV 



is uF u 



10" 



erg cm" 2 s" 1 (see Fig. 11 of MC03a). 



The corresponding absorbed (Ah = 1-94 x 10 cm , 
see MC03a) bolometric (0.02-200 keV) flux is then F pl = 
4.3 x 10 -14 erg cm" 2 s _1 . 



2.4 GX 339-4: ASCA/GIS data 

For the three AS CA/GIS observation of GX 339-4 (see 
IWilms et ail Il999l for a previous study), we also explored 
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the influence of the Nu value on the disc parameters. Again, 
for consistency we fitted the data b y a simple absorbed 
diskbb+powerlaw (|Wilms et al 1 ll999l used a broken pow- 
erlaw, a choice that might influence the radius estimate). 
We chose to probe this effect using three different ways. 
Fir stly by fixing the a bsorption value to the one adopted 
bv lWilms et~aiT(ll999l), i .e. N H , 1 = 0.62 x 10 22 cm -2 (see 
also Zdziarski et al.l 19981 ) , then by fixing it to the weighted 



mean average value we obtained in the following study with 
Swift data (see Section [3]), i.e. AT H , 2 = 0.43 x 10 22 cm~ 2 , 
and finally by leaving it free to vary. Results of those fits are 
given in Tab. [2] 

As noticed with the Swift data, the value of the disc nor- 
malisation tends to decrease whenever the absorption does. 
For example in the third observation the corresponding in- 
ternal radius would reach a value of ~ 17 R g when the value 
of the absorption is fixed to Nn = 0.62 x 10 22 cm~ 2 , and 
drops to ~ 1 R g when it is left free to vary. The disc is not 
even necessary in order to fit the spectra properly when us- 
ing a floating Ah in observation one and two. However, it is 
unclear whether this is due to a real trend or a systematical 
effect linked to ASCA/G1S data analysis. 

A good way of disentangling this problem would consist 
in applying the same processes we did for Swift data (see 
Section [3| . However, as the the number of observations is 
quite low in this case (3 observations compared to more than 
20), we would not be able to draw definite conclusions. As a 
result of this, we chose to keep the results coming from the 
case when Ah is fixed to 0.43 x 10 22 cm" 2 (average value 
obtained with Swift) for the following study. 



to particularly low \ values. We obtain an average value 
of the reduced \ 2 °f about ~ 0.6 on 13 observations. We 
therefore investigated the effect of changing the value of the 
systematics on the calculation of the errors on the fit pa- 
rameters. Changing the systematics from 2.5 per cent to 5 
per cent, the relative errors on e.g. the absorption value in 
obs. 00030009005 jumps from only 1.1 per cent to 1.7 per 
cent. The disc normalisation changes from 3.1 per cent to 4.6 
per cent. Therefore, if the systematics seem overestimated in 
LrPD mode, the actual value has little eff ect on our study. 
We th erefore kept the 5 per cent given in ICusumano et all 
(120051 s ) . 

For certain observations where the disk is dominating, 
the addition of a powerlaw gives an unrealistic photon in- 
dex value (sometimes negative). For those latter cases, we 
did not take into acco unt the 7.5-10 keV energy range (as 
iBrocksopp et al] 120051 do) . The addition of a powerlaw is 
thus sometimes not required. Spectra were also binned in 
order to get at least 20 counts/channel. 

For any given observation of the source, there can exist 
gaps due to Good- Time-Interval (GTI). In the early version 
of the XRT data processing software somewhere the spectra 
for each GTI had to be extracted separately, this effect is 
now taken into account. However, we decided to keep this 
process when fitting spectra, i.e. for one observation we can 
get several spectra and thus different fit parameter values. 
As we wanted to probe any effect on the disc even when quite 
dim, this way of processing allows to check if any change in 
the disc parameters would be due to a temporary instru- 
mental bias. 



2.5 XTE J1817-330, SWIFT J1753. 5-0127, GRO 
J1655-40 and GX 339-4: Swift data 

For XTE J1817-330, we processed the data in the same way 
as in R07 (especially concerning the region sizes for the pile- 
up correction) except that we used the XRT pipeline soft- 
ware version 0.11.5 (2007-08-23 release date) and more re- 
cent response files coming from the CALDB 20071101 ver- 
sion (v009) as well. Of particular importance are the bet- 
ter corrections to the residuals under 0.6 keV in WT mode 
compared to earlier response files. It allows us to perform 
spectral fitting within the 0.3-10 keV energy ra nge though 
adding 3% of systematics (|Campana et al.ll2007h . 

For GX 339-4, GRO J1655-40 and SWIFT J1753.5-0127 
especially in WT mode, the pile-up cor rection was d o ne us - 
ing the second method as described in iMineo et alj (|2007T ) 
and recapped here. For a given observation, different annuli 
(PC mode) or boxes (WT mode) of decreasing sizes were 
extracted in the images and spectra were generated and fit- 
ted with the most convenient model (either a powerlaw or 
a powerlaw+disc). Then when the fit parameter reach the 
asymptotic values obtained with the smallest annuli/boxefl, 
it determines the maximal size of the admitted region. 

For spectra obtained in L rPD mode (GRO J1655-40, 
see also IBrocksopp et al.ll2005T ). spectra have been fitted in 
the 0.5-10 keV range. However, we no t iced t hat the value 
of 5 per cent given bv lCusumano et alj (|2005T 1 tends to lead 



2.6 Fits and flux estimate methods 

Unless explicitly mentioned, all spectra in this study were 
fitted using xspec vl2.3.1ao (|Arnaudl l l996i ). We also used 
the f luxerror tcl script based on Monte-Carlo method that 
was provided by K. Arnaud. It allows the estimation of flux 
errors on each component of the model separately. However, 
we slightly modified the script in order to be able to compute 
the errors within an extended energy range instead of the 
default one. In order to estimate the bolometric luminosity, 
absorbed flux computation has thus been performed between 
0.05 and 200 keV for the hard component. This choice was 
motivated by the fact that: 

(i) the cut-off observed in hard states generally occurs 
around a few hundred keV or the flux contribution over 
200 keV is low for typical soft state (with T ~ 2.5, 

-Fpl, 200.-20000 kev/^pl, 2.-200 keV = 0.1). 

(ii) there should also exist a "cut-off" of the powerlaw at 
low energy (taken into account when fitting by e.g., thermal 
comptonisation models) that can be roughly mimicked by 
the absorption. 



2 In our case, if the photon index differences are less than 5 per 
cent of the asymptotic value. 



A summary of the fit results used in this study for the disc 
geometry analysis (i.e. Sections [4] [5] and [6j are displayed in 
Tabs. \MM and [C] 
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Table 2. Results of the fits by an absorbed diskbb+powerlaw of GX 339-4 ASCA/GIS spectra when changing the absorption value. 





TVh value 


^disc 






x 2 /f 


Obs #1 


0.62 
0.43 

0.32 ± 0.02 (free) 


8.2f 6 X 10 3 
l.lg-ig x 10 3 


0.19 ±0.02 2.3f \ x 10- 10 
0.2lg;27 5.lg :3 x lO" 11 
powerlaw only sufficient 


9.0 ±0.3 x 10-1° 
1.003j;gj$ x 10- 9 


1054.03/954 
1030.70/954 
1025.03/954 


Obs #2 


0.62 
0.43 

0.29 ± 0.01 (free) 


2.4f ° x 10 4 
7.8 40 g x 10 3 


0.18 ±0.01 5.5|-!j! x 10- 10 
0.17 ±0.04 1.4j : * x 10" 10 
powerlaw only sufficient 


1.55 ±0.05 x 10- 9 

1 - 75 i:fo x 10-9 


1279.46/1174 
1201.70/1174 
1184.69/1175 


Obs #3 


0.62 
0.43 

0.27 ± 0.01 (free) 


3.37| g| x 10 4 
2.68 3 ;« x 10 3 

ioc;190 
roo 100 


0.2110J" 1.38i;« x 10-9 
0.279g-pg 3.42| If x lO" 10 
0-46S;'i 1.3lii 4 l x 10-1° 


4.00f° 4 x 10" 9 
4.63 4: g x 10- 9 

5 - 37 i 11 x 10-9 


2970.47/1469 
2506.09/1469 
2366.47/1468 



3 THE HYDROGEN COLUMN DENSITY: TO 
FIX OR NOT TO FIX? 

The value of the hydrogen column density adopted is 
always a central problem when trying to constrain the 
properties of the d i sc. F or many previous studies (R07, 
IGierlinski fc Done. I 120041 ) this value was fixed, either 
by referring to the mos t common adopted values (e.g., 
IGierlinski fc Done. ||2004h . or by taking the one inferred by 
the fit at the highest fluxes. 

Studying the geometrical properties of the optically 
thick disc via X-ray analysis is always a challenge as one has 
to determine the value of two parameters that appear to be 
quite correlated - the column density and the normalisation 
of the emitting disc (that is proportional to its projected 
area). Moreover, this degeneracy will be strengthened by 
intrinsic and/or technical effects; the lower the maximum 
temperature of the disc, the higher the column density or 
the lower the spectral sensitivity of the instrument observ- 
ing, the higher the degeneracy. 

For example, even for sources like GRO J 1655-40 
or GRS 1915±105 that show high disc temperatures 
during their outburst, several authors fix the hydrogen 
column density value w hen using RXTE observations 
(|Gierliriski fc Done. 1120041 for example). This is justified as 
the low efficiency of the PCA under 2 keV would not con- 
strain easily the value of the absorption. However, when 
using observatories as efficient at low energy as Chandra, 
XMM or Swift, there seems to be fewer "a priori" underly- 
ing reasons to fix the Ah for every spectral analysis. 

In this section we will study the effect of letting the 
value of Ah varying versus keeping it fixed, especially dur- 
ing the transition from soft to hard state. We apply this 
to the Swift-XRT data of GX 339-4, XTE J1817-330, GRO 
J1655-40 and SWIFT J1753.5-0127. Initially we evaluated 
the value of Ah at the highest luminosities. In that case, 
its value varies between 4.6 x 10 21 to 2 x 10 21 cm from 
observation number one to 13 in GX 339-4 data, and be- 
tween 2 x 10 21 to less than 1 x 10 21 cm- 2 for XTE J1817- 
330 for example. Then, when a multicolour disc is re- 
quired according to our standards (see hereafter in section 
[4|, we can compute the weighted mean value Ah, wmcan- It 
gives A H , wmcan = 0.43 ± 0.02 x 10 22 cm" 2 for GX 339-4, 
0.12 ± 0.02 x 10 22 cm -2 for XTE J1817-330 (in agreement 
with the value obtained by R07), 0.246 ± 0.011 cm -2 for 



SWIFT J1753.5-0127 and 0.73±0.014 cm" 2 for GRO J1655- 
40. 

If we now plot the dependency between the value of 
Ah and the photon index as shown in the first column of 
Fig. [T] one can see that there seem to be a slight correlation 
between the value of Ah obtained from one observation to 
another. That effects can also be seen when looking at the 
distribution of Ah values presented in the left column of Fig. 

m 

However, as we combine three different models (wabs x 
(diskbb + powerlaw) or wabs x (diskbb + compTT)) with a 
total of five free parameters, one cannot a priori exclude a 
systematic degeneracy among those. Indeed, one can argue 
that the trend between the photon index and the value of Ah 
obtained could be purely due to the model used, as a higher 
photon index would require a high absorption to obtain one 
and the same flux at low energies. 

Therefore, in order to probe if the trend was real or due 
to the fitting process, we proposed to follow the following 
scheme summarised on Fig. [3] Having processed the real 
data as described above (i.e. with a floating Ah), in parallel 
and for each observation we fitted the same spectrum in 
the same way (i.e. adding a disc, or not, when it was not 
appropriate), except that the Ah is now fixed to the 
best weighted mean value obtained above. With those 
new fit parameter values, we then faked a spectrum with the 
XSPEC command f akeit, taking the same response files and 
time elapsed, thus obtaining the so-called fixed then faked 
(ff) spectrum. In a third step, we fitted this ff spectrum 
with a floating Ah, hence obtaining fff (fixed, faked, free) 
parameters. The comparisons in the value of Ah obtained in 
the fff and in the simple case of a floating Ah would allow 
one to determine whether the observed variation is more 
consistent with statistic/systematic deviation or with a real 
trend. 

Results are displayed from Fig. [2] to [S] and Tabs. [3] 2] 
and[5] Two complementary approaches have been used: first 
by comparing both distributions of Ah obtained, and then 
the apparent photon index F-Ah correlation observed. As 
stated in the previous section, both cases where the high 
energy part of the spectrum i s fitted with either a powerlaw 
or a comptt iTitarchuk et all (|l994h model were examined. 
Both analysis seems to be useful, as it seems to lead to 
different conclusions. 
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Figure 1. Column density variation during the transition in GX 339-4, XTE J1817-330, SWIFT J1753. 5-0127 and GRO J1655-40 when 
high energy component in real and faked data are fitted with a powerlaw model. The model is either a wabs*po (o) or a wabs*(diskbb+po) 
(A). Dashed line shows the fit by a powerlaw of the A and the dottcd-dashed line is a fit of all acceptable points (A and o when T > 1). 
See Tab. [3] for the fit results. Note that for some observations a powerlaw component is not necessary to fit the spectra (e.g for GRO 
J1655-40). Therefore the number of points in those graphs can be different from the numbers given in the histograms of Fig. [2] 
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Figure 2. Distribution of TVh when the high energy component in real and faked data arc fitted with a powcrlaw model. Left panels 
correspond to the real data fitted with a floating TVjj, when the right panels correspond to data obtained by faking the real data when 
Nh is fixed to the weighted mean value obtained in the highest states. 
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Figure 3. Schematic process we adopted to probe if the variation in the column density value observed was a systematic effect of data 
analysis or real. This comparison is called the fff method in the text. 



3.1 High energy component fitted by a powerlaw. 

When a powerlaw is used first, one can see in the histograms 
of Fig. [5] the relative differences between the distribution of 
JVh obtained in the real and the faked spectra. As a result, if 
we except the case of GRO J1655-40, the real distributions 
seems to be flatter than the fff ones. We then performed a 
Kolmogorov-Smirnov (KS) test to compare the latter distri- 
butions. The results, shown on Tab. [5] give a low probability 
that the real and the fff spectra could come from the same 
distribution. However, except for SWfFT J1753.5-0127, the 
KS test probabilities obtained are still too high to draw any 
definite conclusions. 

Had the fitted JVh been an artifact of the degeneracy 
with the fit parameter, we would have expected the distribu- 
tion to be similar. To confirm the non degeneracy, we fitted 
the Nh versus F by a powerlaw in both real and fff cases as 
shown on Fig. [T] Note that the choice of a powerlaw is com- 
pletely arbitrary. However, it seems straightforward to probe 
whether there is a systematic trend by examining and com- 
paring the exponent value obtained: a nil value of this latter 
would for example mean that Nh is independent of F. We 
performed this fit using either the whole dataset (dotted- 
dashed curve in Fig. [TJ, or only the observations where a 
disc is required (dashed curve). For the x 2 evaluation when 
fitting, each point in the Ah vs F diagram is weighted ac- 
cording to its error on A^h only. The powerlaw exponent 
values obtained are thus displayed in Tab. [3] 

If we pay attention to the real spectra when a disc is re- 
quired, the exponent values range from a weak dependency 
of 0.15 in GRO J1655-40 to a bigger one of 2.25 in SWIFT 
J1753. 5-0127. The accuracy of these exponent value is higher 
than 2.5a for the four sources studied. If we now examine 
the values obtained with the fff spectra, we note that the 
exponent is still positive, except for GRO J1655-40, suggest- 
ing a possible degeneracy. However, the Ah vs V dependency 
is lower in the fff C£lSG ; ELS shown in in Tab. [3] The exponent 
values obtained in the real case are always higher than those 
obtained in the fff case by a factor 2.5 to 4.7. Examination 
of the same analysis performed on the whole dataset (right 
hand side of Tab. [3} gives similar results. 

We then conclude that if a powerlaw is used as 
model to fit the high energy part of the spectrum, 
the value of Ah seems to vary during the evolution 
of the source, and Sioifi/XRT is able to detect this 



evolution. As a result, we suggest to keep Nh as a 
free parameter when fitting Swift /XRT data with 
this modefj], and, we suggest to use the fff method for 
further Nh variation studies with other observatories. 

3.2 High energy component fitted by 
comptonisation model. 

As the underlying processes resulting in the high energy tail 
are still under debate, we also chose to examine if the pre- 
vious variation of JVh was still evident when using another 
model. We thus chose xspec comptt model, often used as 
a physically motivated alternative compared to a powerlaw. 
But this model has the characteristic to exhibit a low energy 
"cut-off" (Rayleigh- Jeans part of the seed photons coming 
from the black-body emission), which can greatly influence 
our previous analysis. If the "true" emission process is actu- 
ally not extended to low energy, and a powerlaw (which by 
definition exhibits no low energy cutoff) is used, one expects 
the absorption to play the role of an artificial cutoff. Hence 
the Ah value obtained when using a powerlaw will be higher 
than when using comptt. 

Thus we reran the previous analysis with this model 
and the results are displayed in Fig. U and [S] and Tabs. [3] 
and [5] Firstly, we can visually examine fff distributions ob- 
tained on right-hand side of Fig. [4] We note that instead 
of being normally distributed around the adopted weighted 
mean value of Nh, the fff spectra tend here to exhibit a sec- 
ond peak at low A^h for GX 339-4 (the six observations where 
N H ~ 2 x 10 21 cm" 2 ) and SWIFT J1753. 5-0127 (11 obser- 
vations with Nh ~ 5 x 10 20 cm -2 ). This means that faking 
data with a constant Ah can, in certain cases when using the 
comptt model and if we refit the same spectrum with a float- 
ing absorption, lead to quite significant differences in the Nh 
fitted value. Moreover, examining the KS test results (right 
hand side of Tab.[SJ) gives a higher probability that real and 
fff distributions are the same when using comptt compared 
to when using a powerlaw (except for GRO J1655-40). This 
effect is particularly visible for SWIFT J1753. 5-0127 where 
the KS probability reaches 15 per cent. 

In order to probe if any correlation between the state 
and the value of Nh remains, we used the plasma optical 

3 Of course, this will only be valid for high enough Nh- 
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Figure 4. Distribution of JVh when high energy component in real and faked data are fitted with the comptt model. 
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Figure 5. Column density variation during the transition in GX 339-4, XTE J1817-330, SWIFT J1753.5-0127 and GRO J1655-40 when 
high energy component in real and faked data are fitted with the comptt model. See Tab.|4]for the fit results. 
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Table 3. Result of the fit by a powerlaw plotted on Fig.[T] Dashed line of Fig. \T\ correspond to the "only obs. with disc" where dottcd- 
dashed line refers to "all obs." . All uncertainties are 1 <r wide. These results are mainly consistent with a real Ah versus T correlation, 
however partly due to systematic deviation coming from spectral analysis. 







Model : A H 


= axT> 3 




Source 


/3 (only obs. 


with disc) 


P (all obs.) 


real 


fff 


real 


fff 


GX 339-4 


0.47 ±0.05 


0.10 ±0.09 


0.98 ±0.14 


0.09 ±0.06 


XTE J1817-330 


0.85 ±0.32 


0.35 ±0.14 


1.11 ±0.39 


0.14 ±0.10 


SWIFT J1753. 5-0127 


2.25 ±0.38 


0.86 ±0.35 


1.16 ±0.40 


0.15 ±0.09 


GRO J1655-40 


0.15 ±0.06 - 


0.009 ±0.014 


0.23 ±0.11 


-0.01 ± 0.014 



Table 4. Results of the fit by a powerlaw plotted on Fig. [5] Every terms have the same signification as Tab. [3] Here, except with XTE 
J1753. 5-0127, there seem to be less correlation between the value of the warm plasma optical depth and the column density. 



Model : Nn = axr 13 



Source 


j3 (only obs. 
real 


with disc) 
fff 


P (all obs.) 
real fff 


GX 339-4 


-0.02 ±0.01 


-0.025 ± 0.008 


-0.37 ±0.30 


-0.009 ± 0.014 


XTE J1817-330 


-0.003 ± 0.009 


0.008 ± 0.003 


-0.004 ± 0.006 


0.010 ±0.003 


SWIFT J1753.5-0.127 


-0.51 ±0.16 


-0.19 ±0.14 


-0.71 ±0.14 


-1.9 ±0.3 


GRO J 1655-40 


0.072 ±0.031 


0.025 ±0.039 


0.070 ± 0.030 


0.030 ± 0.032 



depth value, r, as an indicator equivalent to the photon in- 
dex (since V depends on the compton parameter y, where 
y = kT e / (m e c 2 ) x max(T, r 2 ) and T e is fixed at 50 keV). Fig. 
[5] illustrates the possible Ah vs t dependency, and fits by 
powerlaws have been performed as in the previous section. 
Values of the exponent obtained for observations where a 
disc is required are consistent with no variation of Ah vs 
r in GX 339-4, XTE J1817-330 and GRO J1655-40. It is 
either due to the fact that the value of the exponent is sim- 
ilar between the real and the fff spectra (GX 339-4 case), 
or because its la error for real spectra is as large as the 
value obtained (XTE J1817-330 case). GRO J1655-40 seems 
even to show inverse behaviour - the value of the hydrogen 
column density seems indeed to correlate with the value of 
the plasma optical depth. On the contrary, SWIFT J1753.5- 
0127 still seems to exhibit significant decrease of Nu when 
t increase. 

Considering the complete dataset (dot-dashed curve in 
Fig. [5] and right hand side of Tab.[4| returns similar results, 
except for SWIFT J1753. 5-0127, where the bimodality of 
fff distribution is quite noticeable, forming two separated 
cluster of points (lower right graph of Fig. [4]). 

Thus, when using comptt to fit the high energy tail of 
the spectra in Swift /XRT data, there is less evidence that 
the absorption column density is varying. As a result one 
should view the fitted value of Ah with caution. However, 
as the statistical errors obtained on r in this study 
are quite large, it is difficult to draw definite conclu- 
sions about the Ah vs t dependency. The simultane- 
ous use of an instrument which gives more accurate spectra 
at higher energy together with our prescribed fff method 
could help to probe any Ah vs r variation. As a general 
conclusion, we suggest to leave the absorption free 
to vary when fitting Swift /XRT spectra. 



Table 5. Results given by a KS test probing the probability that 
the underlying distributions of Ah obtained in Fig.[2]are one and 
the same for the real and the faked spectra. 



powerlaw comptt 



KS test: 


value 


prob. 


value 


prob. 


GX 339-4 


0.34 


0.048 


0.375 


0.051 


XTE J1817-330 


0.41 


0.016 


0.33 


0.09 


SWIFT J1753.5-0.127 


0.73 


7.1 x 10~ 9 


0.26 


0.15 


GRO J 1655-40 


0.28 


0.24 


0.41 


0.017 



4 A RECESSING DISC WHEN Lbol < 10" 2 L E dd 
IN XTE J1817-330? 

Although we have just investigated in which cases one has 
to let the value of Nn float during the fitting process, we 
will examine in this section what is its impact on the disc 
properties is, especially its geometry. Thus, we will still con- 
sider on the one hand, the case when Ah is free to vary, and 
on the other hand, when it is frozen to the value obtained 
at high luminosity (see Sect. [3] for those values). This work 
focuses on the Swift-XRT observations of XTE J1817-330 
which appear to show significant changes in the fitted disc 
properties during its decline. 

Fixing the value of Ah, or leaving it free during the fit- 
ting process has a major impact on the possible detection 
of a disc component. However, by comparing the effect of 
adding a diskbb model, it seems that in both cases the fit 
is improved at high luminosity. In the case of an Nn value 
fixed to the one obtained at high luminosities (as opposed 
to the case where A^h is free to vary), taking the example 
of observation number 13, the \ 2 drops from 13720 to 386.2 



Disc radii variation in black-hole binaries 13 



(respectively 2281.9 to 380.9) when adding the thermal com- 
ponent. 

We now focus particularly on the last four observations 
of XTE J1817-330 and especially observation number 20. 
The goodness of fit summary is shown in Tab. [6] Taking into 
account the degrees of freedom, the situation is less clear, 
it can show that a simple absorbed powerlaw already gives 
a relative good fit when the Nh is free (reduced x 2 = 1-08) 
and adding a diskbb in this case suggests an overfitting of 
the data (the F-test probability is higher than 3 per cent). 
On the contrary, when the Nh is fixed to the value obtained 
at high luminosity, as the absorption in this latter case is 
higher, adding an extra component at low energies, like a 
diskbb, will be required (see Tab. [6}. This effect is due to 
the fact that, as we demonstrated in previous sections, when 
using a powerlaw for the high energy component, the ab- 
sorption seems to slightly decrease when entering the harder 
states. Thus an extra component at low energy such as a disc 
would not be required anymore. On the contrary, if we fix 
the Nh to a high value, then a disc component would be 
necessary to compensate the flux lost by the absorption. 

In the next sections, considering that this effect could 
affect the possible detection of a disc component, we adopted 
the following strategy when analysing Swift data: 

• For easier comparisons and consistency with earlier re- 
sults given in the bibliography, we fitted the high energy 
part of the spectra with a powerlaw. 

• If ever the spectrum is well fitted by a simple absorbed 
powerlaw (threshold chosen: \ 2 ' l v < 1-2), or if adding a disc 
component to the latter model gives a high F-test probabil- 
ity (hereafter, > 10 -2 ), then no disc component is added. 

• On the contrary, when an extra component is necessary, 
we study the evolution of the parameters and we allowed the 
absorption to vary (unless it is explicitly mentioned as for 
XTE J1817-330 case where both varying and fixed Nh case 
were studied). 

We can also compare our results with those obtained 
by R07. For example, in observation number 19 we obtain 
X 2 jv = 18.33/23 with a simple absorbed powerlaw only. 
The absorption value obtained in this case is quite low 
(Nh = 7.2 ±5 x 10 20 cm" 2 ), however con sistent with earlier 
Chandra observation dMiller et al ] |2006ch . Even if we set the 
iV H value to 1.2 x 10 21 cm" 2 , we obtain /v = 20.02/24, 
which is already sufficient to model the spectrum properly 
(the addition of a diskbb does not improve the fit). It is not 
clear why R07 added an extra component in that observa- 
tion. We have similar concerns for obs. 20 and 22. Examining 
the x 2 values obtained (see Tab. [SJ , one can see that adding 
a disc component is not necessary for the four last observa- 
tions except, as noted above, in observation number 20. Note 
also that in obs. 20 the spectrum is only made of 27 chan- 
nel bins. Performing a F-test between both models (with or 
without disc) gives a probability of 2 x 10~ 3 , which is still 
quite high and still cast doubt on the reliability of adding 
such thermal component. This disc measurement was how- 
ever kept for this section's study. 

The variation of the hydrogen column density value has 
also been reported in previous papers involv ing observa- 
tions with Swift (see e.g. iBrocksopp et al]|2005l for the GRO 



J1655-40 outburslQ) for objects in their soft state but also 
in their hard states (jOosterbroek et ahlll996l in V404-Cyg). 
We also saw in Sect. [5] that ASCA/GIS observations of GX 
339-4 can be interpreted in the same framework and lead to 
different conclusions concerning the absence or the presence 
of a thermal component whether the value of Nh is fixed or 
not. 

We now focus on the remaining observations of XTE 
J1817-330 at higher luminosity, (observations number two 
to 17 and even 20 if Nh is fixed). When a disc is required, 
one can track the evolution of the source as its flux is 
monotonicall y declining from the first to th e last observa- 
tion (see also Gicrliiisk i. Done fc Pagell2008D . As noticed in 
iGierlinski. Done fc Page! ( 20081 ) by contemporaneous obser- 
vations with RXTE, in the hardness-intensity diagram the 
source is fading from the high soft state (HS) to the low 
hard state (LH) by tra nsiting through the lowe r intermedi- 
ate state branch (IM) (|Homan fc Bellonil 120051 ) . Note that 
observations number 15, 16 and 17, of particular interest as 
we will see, are located in this latter branch. 

The whole detailed evolution of the disc properties is 
displayed in Fig. [5]when using a simple powerlaw as a model 
to fit the hard component of the spectra, and in Fig. Qwhen 
using comptt for the same purpose. As mentioned before, 
we also examined both the cases where Nh is fixed or not: 
as there are 4 observables (Ldisc, Lboi, Rin, kT[ n ), there are 
in theory six 2-dimensional plot combinations. However, for 
simplicity, we restricted this number to the 4 most important 
for probing the source evolution and we discarded the Lboi 
vs Ldisc and the kT ln vs Lboi plots. 

Examination of the first row of graphs in Fig. |6](fcTi n vs 
Ldisc) seems to clearly show an obvious trend: where at high 
luminosity the temperature seems to decrease monotonically 
as a unique powerlaw, the last three "double points" (there 
was two GTIs for each of those three observations) shows a 
clear drop in temperature. We then tried to fit this kT- ln vs 
Ldisc relationship with two models: first using a simple pow- 
erlaw, and second with a broken powerlaw. The f it algorithm 
uses t he standard merit function described in IPress et alj 
(|1992| ). § 6 .7. It does take into account both the errors 
on the Y (inner temperature) and X (luminosity) axis: the 
principle of this fi tting proceedure is described in details in 
IPress et all (|l992l ). § 15.3, "Straight-Line Data with Errors 
in Both Coordinates". This algorithm was also used when 
fitting the kT ln vs -Rin and Lboi vs i?; n relationships. 

The results are given in the first set of rows of Tab. [7] 
first, a simple powerlaw to model those data give a very bad 
fit, whenever the value of Nh is fixed or not (Xspi/^ ~ 5-0 or 
~ 10.) . Using a broken powerlaw significantly improves the 
fit as it drops to a value of 1.15 (Nh free) or 2.3 (Nh fixed). 
Second, the upper flux powerlaw index (0.27 or 0.26L0.02) 
is very close to the value expected for a disc emitting with 
a constant inner radius (fcLl n oc L°;2). By comparison, R07 
obtained a powerlaw index equal to 0.233 ± 0.006. It is in- 
teresting to note that the second powerlaw index found is 
quite steep and close to 0.77 (Nh free) or 0.6 (Nh fixed) at 
lower flux. The errors on those values, although bigger, are 
still consistent with a break which would occur at around 
0.014 Lsdd (with orbital parameters and distance given in 



Based on the same data set used in this paper. 
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Tab.[T}. As the temperature seems to suddenly decrease for 
a rather constant disc luminosity, one expects its normali- 
sation and hence the internal radius to increase in order to 
compensate for this fall. This is what is examined in the six 
lower panels of Fig. [6] In observations 16 and 17, the inter- 
nal radius values double and triple from those measured at 
higher luminosities. 

The second row of Fig. [6] shows the variations of fcTm as 
a function of Ri U . A broken powerlaw has been used again 
to fit the data. Since for a pure multicolour disc each ring is 
dissipating the gravitational energy available in radiation as 
it is optically thick, it can be used as a direct probe of the 
accretion rate. In that case we have: 



Tin = 



/ 3GMM d 



1/4 



(•>) 



\S-KR? n a B 

(see e.g lBelloni et al.lll997h . where Md is the mass accretion 
rate flowing through the disc only. This supposes that 100 
per cent of the mass coming from the donor star is accreted 
entirely through the disc until Ri n . Therefore, renormalis- 
ing the same formula to Eddington accretion rate and using 
I/Edd = ??EddMEddC 2 , we obtain: 



1 keV 



= 5.31 



M 
M© 



-1/4 



m d 

»?Edd 



1/4 



Rn 



-3/4 



if) ' ^ 



where riid = Md/Afodd- Note t hat rhd ^ ihtotai in the 
general case (see for example Ivan der Klisll200l] for 
evidences of this discrepancy in LMXB). 

Hence, parallel lines with a slope of (-3/4) can be drawn 
in a log-log temperature versus radius diagram and the nor- 
malisation value only depends on the values of the mass, 
the radiative efficiency at Eddington luminosity and the ac- 
cretion rate. This has been done in the second rows of Fig. 
[6] and [7] Below the apparent break in the relationship, we 
observe that the disc seems to decrease its accretion effi- 
ciency by increasing its radius, keeping its accretion rate 
constant, especially when Ah is kept free. This is further 
supported by the value of the slope given by the fit even if 
the errors on it are quite high (—0.70 ±0.33 consistent with 
-3/4). By taking 7/Edd = 0.1, we conclude then that the ac- 
cretion rate is close to 0.005 Eddington (±0.0025 with the 
large uncertainties on the BH mass). This effect also seems 
to hold when using a comptt (see second row of Fig. [TJl . Us- 
ing Ldisc = 27r-Rf n o"BT]' n , we can easily calculate the value of 
the accretion rate through the disc0, i.e.: 



m d 



2 / Ldis 

o^Edd 



Rin 

R E 



(7) 



Moreover, as variations in the mass ratios of BHBs are 
far less than those on accretion rates (10 s ), one can set a 
typical mass value (e.g 8Mq) and Eddington accretion effi- 
ciency (e.g r^Edd = 0.1) in order to estimate the evolution of 



5 Note that the equation number [7] is obtained in the framework 
of the multicolour disc model. One known limitation is that it 
gives a quite high accretion efficiency when the radius approach 
the ISCO, as in that case r](R in ) diskbb = 3/2(R g /_R in ) = 0.25 ^ 
0.06 when _R; n = 6 R g for example. 



the accretion rate among different sources. This will be the 
subject of Section [5] 

The beginning of the disc recession can also be noticed 
in the two last rows of graphs of Fig. [5] and [7] and in Tab. [7] 
for the results of the fit by a broken powerlaw: the crucial 
transition observed in the geometry for this source seems 
then to occur at the same bolometric luminosity level (~ 
0.015I/Edd)- This is due to the fact that until observation 
number 15, the disc contributes to more than 90 per cent of 
the overall luminosity (it drops to 58 per cent in observation 
17). 

We also compared and included with our results the 
XMM observation of the source analysed by (|Sala et al.l 
120071 . S07). For a better comparison with our data process- 
ing manner, we used the fit results that S07 obtained with an 
absorbed powerlaw+diskbb. It appears that the disc normal- 
isation value is still in the range we obtain with Swift data 
in the highest fluxes state (as -ffdisc — 2000) although the 
value of the photon index differs quite a lot when the value 
of TVh is tied to a value of 1.5 x 10 21 cm -2 . With T ~ 1.6, the 
XMM-Newton observation seems to exhibit a non-recessed 
disc in a low-hard state, similarly to the observations anal- 
ysed by M06 on GX 339-4. However, if we focus on the value 
of the fluxes obtained using the fit parameters values given 
in S07 with Ah frozen, it results in an unabsorbed disc flux 
of fdisc — 1.1 x 10~ 8 erg cm -2 s~ x and a high absorbed 
powerlaw flux of F pl = 2.4 x 10~ 7 erg cm -2 s -1 in the 0.02- 
200 keV range. Thus in terms of Eddington luminosity it 
gives Ldisc ~ 0.060 L E dd but L bo i ~ 1.3 L E dd (M = 6M , 
d — 6.3 kpc, see Tab. [JJ. Unless one cannot rule out the 
possibility that the source lies in our vicinity (i.e d ~ 1 kpc, 
in that case Lhard ~ 0.012 LEdd); the overall luminosity is 
quite high. 

We can also note that S07 obtained a much higher pho- 
ton index (F ~ 2.3) when letting the value of Ah free to 
vary, but the high energy flux is still completely dominat- 
ing the spectrum. Such a high bolometric flux is quite puz- 
zling as the closest Swift /XRT observations of the source 
taken 6 days before and two days after (7th and 15th 
March 2006, corresponding to ObsIDs number 03 and 04) 
are one order of magnitude less luminous (Lboi ~ 0.13 and 
O.llLEdd respectively). We note however that using comptt 
to model the high energy part of the spectrum gives a 
more consistent bolometric luminosity value, as in that case 
Lboi ~ 0.06LEdd- The discrepancy could come from the fact 
that the x 2 value is far better with this latter model rather 
than when using a simple powerlaw {x* /v = 1.19 instead 
of 1.45). However, in spite of those drawbacks, for the con- 
sistency and comparisons among the different sources and 
observation, we kept these results given by the fit with a 
powerlaw (with .Ah free). 

In the case that XTE J1817-330 lies at lkpc, we would 
also obtain R ln ~ 0.77±0.05 R g < Rs in the highest state of 
SWIFT data (as estimated from the first five observations) . 
This casts doubt on such a low value for the distance. In 
order to get at least R[ n ~ Rs, one has to set a distance 
of 2.6 kpc and we thus took this model dependant value as 
minimal distance to XTE J1817-330 for the following multi- 
source study. Moreover, if we choose d — lkpc, the beginning 
of the disc recession observed in SWIFT data would then 
occur at about 4 x 10 -4 LEdd- 
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Figure 6. Variation of the disc properties during the decline of XTE J1817-330 when fitting the hard component with a powerlaw. The 
A;Ti n vs I/disc/^Edd, kT[ n vs -Rin/Rg and ^bol/^Edd vs fiin/Rg data sets have been fitted by a broken powerlaw (solid lines, see text 
and Tab. [7] for the parameters of the fit). The dotted dashed lines plotted in the second row of graphs are showing constant rhd profile. 
From upper right to lower left: rh^ = 0.002, 0.005, 0.01 with rygdd = 0.1. 
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Figure 7. Variation of the disc properties during the decline of XTE J1817-330 when fitting the hard component with the comptt model. 
Here again the beginning of a disc recession seems to be needed. Even if the behaviour at high luminosity is less constant. Same values 
as in Fig. [6] for m d in the second graph row. 
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Table 6. x 2 values when fitting the 4 last Swift observations of XTE J1817-330 with or without a disc component and using the latest 
response files. Even when Nn is fixed to 0.12 X 10 22 cm -2 , the reduced x 2 is quite good enough (i.e. < 1.15) to model the data with a 
simple powerlaw in observations 19, 21 and 22. Adding a disc component would overfit t he data in this case. N ote that the values of the 
f-test probability are just informative as its use is controversial for such small data sets llProtassov et al.ll2002l) . 







x 2 






Model 


Obs. 19 


Obs. 20 


Obs. 21 


Obs. 22 


wabs X po 


±|f^ = 0.80 


2|f = 1.08 


2Mi = 0.90 


2§|° = 0.93 


wabs X (diskbb+po) 


= 0.83 


TS 8 = °- 87 


aya = o.9o 


= 0.80 


ftest probability 


0.63 


0.036 


0.39 


0.022 


wabsxpo (A'h fixed) 


3%g2 = 0.83 


3M1 = 1.54 


36.82 _ i i c, 


39-81 - I 9 
39 ~~ LMZ 


wabs X (diskbb+po) (JV H fixed) 


iff^ = 0.83 


^f- = 0.97 


= 0.93 


= 0.89 


ftest probability 


0.37 


0.002 


0.015 


0.038 



Table 7. Fit results of Fig. [6] 





kT in = 
kT in = 


N x (L diBC /(L Edd x 6)) ai when L disc /L Ed d 
N x (L d isc/(iEdd x b)) a 2 when L disc /L Edd 


> b 
< b 


kT in = N X (L diBC /L Edd ) a ^ 




b 


N ai 02 


Xbknpl/^ 




Nft free 
Nn fixed 


0.014 ±0.001 
0.014 ±0.001 


0.46 ±0.021 0.27 ±0.018 0.77 ±0.15 
0.47 ±0.016 0.26 ±0.018 0.57 ±0.1 


22.9/20 
46.04/20 


109.2/22 
214/22 




-Rin/Rg 


= N x (L bol /(L Edd x b))"i when L bol /L Edd 
= N x (L bo i/(L E dd x b))"2 w hen L bol /L Edd 


> b 
< b 


Rin/Rg = N x (L bol /L E dd) a3 P' 




b 


N a± ct2 




x 2 p i/^ 


Nn free 
Nh fixed 


0.015 ± 0.003 
0.015 ± 0.002 


5.56 ±0.32 -0.05 ±0.03 -2.4 ±1.4 
5.20 ±0.30 -0.02 ±0.02 -1.72 ±0.6 


55.8/20 
37.09/20 


152/22 
299/22 






Ki n /R g = N x (kT in /b) Q i when kT in > b 
R in /R s = N X (kT in /b) Q 2 when kT in < b 




R in /R g = N x (kT in ) Q =pi 




b 


N ai a2 


XbknpA 




Nn f ree 
Nn fixed 


0.44 ±0.10 
0.50 ±0.03 


5.20 ±0.60 -0.12 ±0.15 -1.44 ±0.7 
4.96 ±0.57 -0.02 ±0.09 -0.99 ± 0.2 


82.7/20 
47.7/20 


107.34/22 
175.6/22 



5 MULTI-SOURCES COMPARISONS 

While XTE J1817-330 is one of the best-studied examples, 
we need to see how typical its behaviour is for BHBs. We first 
extended the previous study to the other available SWIFT 
data of candidate black holes, i.e GX 339-4, SWIFT J1753.5- 
0127 and GRO J1655-40. According to the results obtained 
in Sect. [3J we let the value of A^h free to vary and used an 
absorbed powerlaw (and a diskbb when needed) to model 
the spectra. Observations where a thermal component is re- 
quired are plotted in Fig. [5] 

Contrary to the case of XTE J1817-330, there is less sys- 
tematic evolution of the disk radii. Indeed, if a faint increas- 
ing of its value seems to occur at low luminosity in GX 339- 
4, the statistics are not sufficient to draw any conclusions. 
The additio n of the broadb a nd an alysis of lower luminosity 
states led bv lTomsick et al.l (|2008l ) (T08) changes this a bit. 
Indeed, when using reflection models, T08 seem to obtain a 



geometry that does not evolve through the hard state (the 
solid angle value remains constant and high enough). On the 
contrary, the normalisation of the multicolour disc tends to 
exhibit quite high value, whatever the high energy model 
chosen. For purpose of simpler comparison between sources, 
we included the results coming from their modelling by a 
simple absorbed diskbb+powerlaw. Moreover, the evolution 
of the bolometric luminosity is small (changes by a factor of 
~ 5 only) and thus further observations of the decline would 
be necessary. GRO J1655-40 exhibits unusual behaviour; al- 
though the fitted radii change by a factor of two, the rela- 
tionship between the temperature and the disc luminosity is 
not monotonic but rather quite complex. SWIFT J1753.5- 
0127 even shows a disc that tends to get hotter when its 
luminosity decreases, with a higher internal radius in more 
luminous state. However, we were limited by the statistics 
and narrow range in luminosity sampled. 
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GX339-4, M = 7. Mo, d = 8kpc 



GRO J1 655-40, M = 6.3 M , d = 3.2kpc 



Swift J1 753.5-0127, M = 6 M , d = 2.9kpc 
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Figure 8. Variation of the disc properties of GX 339-4 (first column), GRO Jf655-40 (second column) and SWIFT J1753.5-0127 (third 
column) as recorded by SWIFT. Note that in the case of GRO J1655-40, the source is in outburst whereas it is mainly in the decline 
for both others. On top left panel, fit of the cloud of point by a powerlaw gives an index of ~ 0.10. For SWIFT J1753. 5-0127, this index 
is found to be negative (~ —0.21): the disc luminosity decrease when the temperature increase. Concerning the fcT; n vs ij; n plots, same 
coding in ihj as in Fig. [6] has been used. For GX 339-4, both points corresponding to the observations analysed by T08 (▼ symbols) with 
a simple absorbed diskbb+powerlaw have been added. 



We then compared the whole set of sources in a single 
study. This has been done by overplotting on the same graph 
the previous results, and adding the few points obtained at 
lower luminosities. These latter points mainly come from the 
low efficient BHB XTE J1118+480. However, other recent 
or older observations of the previous sources (T08 for GX 



339-4 and S07 for XTE J1817-330 as mentioned before) have 
been added. The overall results are displayed in Figs. [9l 1101 

nu nana in] and mi 
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Mass, distance and inclination angle uncertainties 

The uncertainties plotted directly on the graphs are the 
statistical errors arising from the fits. We also attempted 
to estimate the systematic uncertainties due to the the or- 
bital parameters (see Tab. [JJ. For simplification, we com- 
puted the lower and upper boundaries by taking the limit 
values in masses, distances and angles. For example, as 
the disc luminosity goes with I/disc/iEdd oc Fdi B cd 2 /(M x 
cos(i)), all the disc luminosities will be shifted by a factor 
dmax/(-^min x cos(i max )) if one wants to obtain the highest 
value of Ldisc / Lsdd • 

This corresponds to a vector in Fig. [9] We note, how- 
ever, that this method is somewhat overestimating the shifts 
in the corresponding diagrams as there may be strong cor- 
relation among those parameters. As the mass is often de- 
termined via the mass function, in that case M x sin 3 (i) 
must remain constant when evaluating the uncertainties. In 
addition, since for the recently discovered black-hole candi- 
dates (J1817-330 and J1753. 5-0127) there is no proper or- 
bital measures yet, we have set quite large uncertainties. 



5.1 Disc geometry evolution 

We first investigate the global relation between the disc lu- 
minosity and its inner temperature (Fig.|5J). As was pointed 
out in the previous sections, a constant truncation radius 
would mean a constant relationship with Tin oc L^isc- We 
thus have drawn in Fig. [5] this relation, normalised to the 
temperatures obtained in the highest luminosity states of 
GX 339-4 and XTE J1817-330. An examination of the dia- 
gram shows that around Ldi sc ~ 10 _2 Z/Edd, the points seem 
to slightly deviate from this relationship toward the lower lu- 
minosities. If we then include the points at lower luminosities 
and in quiescence (mostly coming from XTE J1118+480), 
this deviation appears more clearly. The fact that the points 
are located below the relationship would thus indicate an in- 
crease of the multicolour disc normalisation iodise, hence the 
radius (normalised to R g ). 

As pointed out in the previous paragraph, GRO J1655- 
40 (• points) exhibits a rather different behaviour as for 
the same disc (and even bolometric, see Fig. I10|l luminos- 
ity, the inner temperature is higher in the highest states. 
This is more obvious when considering that their points are 
not located around the dashed relation in Fig. [5] The bias 
coming from the distance and the mass uncertainties could 
explain part of this discrepancy. However, as the disc lumi- 
nosity should be about 10 times higher in order to fit this 
relation, it implies changes of a factor three in distances or a 
factor 10 in masses between sources (or some combination of 
the two). Another reasonable explanation could be that the 
truncation radius value in the highest states of GRO J1655- 
40 is closer to the horizon than in the other X-ray binaries 
studied, increasing its efficiency for an equivalent accretion 
rate. This may imply higher spin. 

The apparent overall increase in radius when the lumi- 
nosity decreases was then probed by the relation between 
the disc or bolometric luminosity (Figs. [TOl and [Til and the 
inner temperature (Fig. [12} or the photon index (Fig. I13[) . 
When examining the disc luminosity as a function of the in- 
ner radius (Fig. II 1 ft . this trend seems to be confirmed with 
a few exceptions. Among them are some of the GX 339-4 
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observations, and obviously the one observation from XMM 
by M06 (A point). 

As pointed out by these authors, the radius value in- 
ferred by the thermal component is close to 2 R g , even 
though the disc luminosity is very low (Ldisc ~ 10 _3 I/Edd), 
suggesting a non recessed disc even in the low accretion 
rate states. However, considering the bolometric luminos- 
ity (see e.g Fig. I10[) , we note a different behaviour. In that 
case we have Lboi ~ 1.4 x 10" I/Edd- This suggests that the 
source was in quite a high state during this observation, even 
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code is the same as in Fig. [9] 



though being spectrally hard (as the photon index F ~ 1.46). 
Indeed, as shown in Fig. 1101 and looking at whole data set, 
the evolution of Lboi vs R[ n seems to be less discrepant and 
an increasing radius is consistent with Lboi < 10 _2 LEdd- 

This is emphasised by the analysis of Fig. 1121 as the 
inner radius values seem to roughly anti-correlate with the 
inner temperatures when kT[ n < 0.5 keV. The data seems 
to follow one of the dashed lines for the relation between T| n 
and Ri n (Eq. [6]l for a 8M© BH and a constant accretion rate 
in the disc, m^. 

However, by using (Eq. [7} and taking into account the 
different BH masses, we infer the evolution of the disc ac- 
cretion rate rhd- As shown in Fig. 1141 the disc accretion rate 
seems to roughly scale with the bolometric luminosity in 
the highest states, i.e Lboi > 10 _2 LEdd- The behaviour un- 
der this value is less clear. On the one hand, XTE J1817-330 
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Figure 13. Photon index T vs the internal radius of the optically 
thick disc. Colour/symbol code is the same as in Fig. [9] 
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Figure 14. Bolometric luminosity as function of the accretion 
rate in the optically thick disc obtained for our sample, using 
cquation[7] Points arc subject to global horizontal shifts due to the 
actual Eddington efficiency value used. The set of curves shows 
theoretical relations. Dashed : Lbol ~ ^disc °c when i?; n = 
6 R g and the efficiency at Eddington luminosity r?Edd = 0.1. 
Under ~ 10 _2 LEddi continuing relationship when Lbol °c m^. 
Supposing rh tota i ~ m d we have: ADAF solutions (dotted), MH02 
(Dot-dot-dashed) and RC00 (Dot-dashed). The value of a and 
p are the same as in Fig. 1161 Errors on m d were computed by 
propagating those on the temperature and the inner radius only. 
Colour/symbol code is the same as in Fig. [9] 



seems to exhibit a beginning of rapid decay of the radiation 
efficiency (Lboi divided by a factor 2 for a rather constant 
rhd between observation number 15 and 17). On the other 
hand, considering the whole data set at those low luminosi- 
ties, the obvious lack of disc detections and the larger un- 
certainties are insufficient to draw the profile of Lboi as a 
function of rhd with accuracy. It is however compatible with 
Lboi/LEdd oc m 2 by taking into account the sources observed 
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Figure 15. Disc fraction luminosity diagram of the whole set 
of observations presented in this paper. The arrow represent the 
supposed path in the DFLD during an outburst, starting from 
quiescence. Colour/symbol code are the same as in Fig. [5] 

in quiescence. This corresponds to a scenario in which a re- 
ceding, but radiatively efficient, outer disc feeds an inner, 
radiatively inefficient flow. 

However, the canonical BHB states do not de- 
pend on luminosity (see e.g th e classification scheme of 
iMcClintock fc Remillardl l2003bl which is based on spec - 
tral, timing and radio behaviour or iHoman fc Bellonill2005h . 
Therefore we plot the evolution of the inner radius as func- 
tion of the photon index, F in Fig. 1131 The evolution is 
then as follows. In the softest states, radii remain quite 
constant, as for example the average value obtained for 
XTE J1817-330 and GX 339-4 when T > 2.4 (14 obser- 
vations) is -Rm, avg = 4.91 R g , with a standard deviation of 
<jR in rma = 0.59 R g . On the contrary, for lower photon in- 
dices, the data are compatible with an increase of the inner 
radius when the spectrum hardens, as long as quiescence is 
not reached. Indeed, in quiescence the photon index values 
are usually higher (r ~ 2). This effect was noted in other 
sourc es as well (e.g. in V404 Cyg F ~ 2.1, see lNaravan et all 
1 19971 ). 

Finally, we investigate the evolution of the relative 
strength of each component (i.e emission from the opti- 
cally thick vs optically thin medium) during the outburst. 
"Hardness-Intensity Diagrams" (HID) are commonly used 
for such analy sis because of their emission model indepen- 
dence (see e.g iBelloni et al.l [2004) . However, they are not 
easily adapted to multi-instrume nt and multi-object analy- 
sis. Fo llowing the method used bv lKording. Jester fc Fender! 
(|2006t ). we use the so-called "Disc Fraction Luminosity Dia- 
gram" (DFLD) that aims at generalising the concept of HID: 
in HIDs harder states are located on the right of the diagram 
and correspond to higher powerlaw contribution. For better 
ressemblance, in DFLD we use on X axis the "Power law 
fraction" (PLF): 



which describes well the disc contribution as the "Disc Frac- 



tion" (DF) is then in turn: 

DF = 1 — PLF. (9) 

The DFLD corresponding to our data set is plotted on 
Fig. [15] Of primary importance is first the pattern drawn 
by XTE J1817-330 during its decline. As the bolometric flux 
fades, the optically thick component is dominating and even 
seems to increase its contribution (PLF falls by a factor 2 
between observation number three and 13). Then the PLF 
suddenly increases and this turnover corresponds to the ob- 
servations (number 15, 16 and 17) where a recessed disc be- 
gins to appear. However, it is not accompanied by a drastic 
change in bolometric luminosity, suggesting that the inner 
part of the disc has been replaced by a quite efficient corona. 

We also note that the powerlaw fraction value in qui- 
escence is equivalent to that observed in the softest states. 
Therefore we suggest that the track followed by a BH X-ray 
binary should be roughly equivalent to the one plotted on 
Fig. 1151 (black arrow). From quiescence, the source will pro- 
gressively increase its PLF as it usually enters in outburst 
by its low-hard state. Then, the powerlaw fraction decreases 
due to the presence of the thermal component in the for- 
mer "high-soft" state and finally returns to the hard and 
quiescent states. 

However, this simple pattern is susceptible to be altered 
whenever a high intermediate state with a strong powerlaw 
component is present. In order to confirm this trend, further 
monitoring is therefore still necessary. 

5.2 Comparison with timing semi-empirical 
relationship and accretion models 

5.2.1 ADAF and inner disc evaporation models 

Further investigations were then led to compare our results 
to available accretion models. Our goal is to constrain the 
mechanisms responsible for the optically thin and radia- 
tively inefficient flows. F or that purpose, we mainly used the 
mode ls proposed by (|Czernv. Rozariska fc Kuraszkiewicj 
|2004 CRK04) as particular attention is paid to the evo- 
lution of inner disc radii. We compared our data to the fol- 
lowing three model s - the classical strong ADAF solutions 
( Abramowica Il995| = "A" model ), and two other models 
( Meyer fc Mever-Hofmeisterll2002l = "B" model or "MH02" 
and iRozariska fc Czernvl |2000| = "C" model or "RC00") 
where evaporation of the inner part of the cold disc is taken 
into account. In these latter models, the evaporation effi- 
ciency is driven by conduction between the corona and the 
disc in presence of a magnetic field. 

We also identified the evaporation radius R OV ap de- 
scribed in those papers with the internal radius we obtained 
from our fits (R C va P = Rm). As modeled by CRK04, one can 
then infer the dependency between Ri n and Lboi simply by 
tuning the values of two parameters, the classical a viscosity 
parameter and the magnetisation f3 = -P g as/(-Pma g + -Pgas) in 
the hot plasma. 

We then used the equations given in CRK04 and 
briefly explored the (a, f3) plane for models "A" and "C" 
(P plane only for model " B", as the analysis driven by 
iMever fc Mever-Hofmeisteil was only computed with a = 
0.3). One constraint we assess is the values of the luminos- 
ity when the inner radius seems to begin to increase. We set 
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Figure 16. Same as Fig, 1101 but with models overplotted. Con- 
tinuous curve refers to the accretion model or timing relationship 
used (see text). Parameters value used: ADAF (dotted curve) a = 
1.33 X 1CT 2 and /3 = l.e - 8. MH02 (Dot-dashed curve): a = 0.3 
and/3 = 0.11. RC00 (Dot-dot-dashed): a = 3xl0 -5 and j3 = 0.11. 
For the timing relationships, we use either v\ = i^Keplcrian(-Rin) 
(dashed line) or v\ = iAjound(-Rin) with h^/Ri n = 0.1 (solid line). 
Colour/symbol code is the same as in Fig. [9] 



this value arbitrarily to the one obtained with XTE J1817- 
330 (see Tab. 0, i.e L bol ~ 0.015 L Ed d and R hl = 5.6 R g . 
Each model with an example of (a, 0) values that fulfil this 
requirement was then overplotted on Fig. 1161 

For the ADAF ( "A" ) model, the adopted value for the 
viscosity parameter was ct = 0.013. As shown as well in Fig. 
2 of CRK04, the radiative efficiency depends only slightly on 
the value of the magnetisation parameter in this case. De- 
creasing (3 only tends to deviate a little towards the lower 
luminosity the Lboi vs Ri n curve for large values of the inner 
radius. As a result, even with such a low value as (3 = 10~ 8 
(i.e plasma highly magnetised) , the curvature does not al- 
low the ADAF model to reproduce the very low luminosity 
observed in quiescence if the viscosity parameter ct is the 
same between the highest and quiescent state (as shown 
on Fig. I16p . Any attempt to fit the points in quiescence 
would require a value of a ~ 3 x 10~ 4 , however the state 
where the radius would begin to increase would then occur 
at Lboi ~ xl0~ 5 L E dd- We note that the ADAF model fits 
well the beginning of the recession observed in XTE J1817- 
330. This model predicts indeed that 

R\n 4 • -2 /ir,\ 

— — oc a m . (10) 
R g 

Moreover, in the highest states, the accretion process 
still remains quite efficient, giving Lboi/LEdd oc m (see Eq. 
16, 17 and 18 of CRK04 for further details). We therefore 
get: 



Rg \ ^Edd / 

It is interesting to note that this exponent value of - 



2 is consistent with that obtained when fitting by a broken 
powerlaw the Ri n vs Lboi points for XTE J1817-330 (see Tab. 
0: «2 = 2.4 or 1.7). Moreover, when extended to larger radii, 
the ADAF solution seems to be in quite good agreement 
with the data obtained with EUVE for the XTE J1118+480 
outburst as well. 

Nevertheless, we have shown that the ADAF solutions 
do not seem to fit the Lboi vs R[ n relationship in both high 
and quiescent states. On the other hand, the approximate 
relation giving the accretion radiative efficiency (equations 
16, 17 and 18 of CRK04) seems to fit the data, as shown 
on Fig. [14] (dotted curve) and if we assume that rh scales as 
rh d . 

Compared to the ADAF model, the models where evap- 
oration of the inner part of the disc is taken into account 
seem to be a better overall description of the Lboi vs R ln evo- 
lution. In these cases, the magnetisation value has far more 
impact. Taking (3 = 0.11 (i.e P mag ~ 8P gas ) with the MH02 
model allows us to fit the overall shape of both the quies- 
cent and the low state observations (except XTE J1118+480 
EUVE data). A similar result can be obtained with RC00 
model with a similar value of the magnetisation, but quite a 
lower value of the viscosity parameter (a = 3 x 10 -5 ). The 
main drawback concerning such a low value of a can be seen 
on Fig. [14] (dot-dot-dashed curve). The radiative efficiency 
does not drop at low m and the accretion flow remains ra- 
diatively efficient, inconsistent with the points in quiescence. 
No other (a, f3) combination has been found which better 
mimics the behaviour in both Lboi vs Ri n and Lboi vs m 
diagrams. This problem does not seem to occur with MH02 
model (dot-dashed curve) where the value of a — 0.3 is 
higher. 



5.2.2 Using the Variability Plane 

We subsequently tried to compare our spectral fit- 
ting results with the expected timing behaviour. Indeed, 
iMc Hardy et all (|2006l ) found that AGN and stellar soft 
state BH could be unified by taking the accretion rate and 
the BH mass into account. Indeed, a fundamental plane 
seems to link the typical frequency observed in the power 
spectra (especially to the lower high frequency Lorentzian 
v\), the BH mass (M) and the total accretion rate rh 




iKording et all (|2007h found afterwhile that this rela- 
tionship could be extended to hard states, even if the nor- 
malisation is not exactly the same, neither is the character- 
istic frequency (in soft state of BHB and AGN, it is the high 
frequency break). If typical frequency observed is linked to 
a Keplerian motion near the internal radius of the accretion 
disc Ri n , we get: 

ia cx (M) 1/2 (i? in r 3/2 cx (M)- 1 (^) ^ , (13) 

when normalising Ri n to the gravitational radius R g oc M. 
Thus replacing v\ in equation[l2]by the previous value gives: 
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-3/2 



(M) oc rh 



Ri 



-3/2 



For inefficient accretion flows that may occur at low 
luminosities, one can also have: 



Lhol . 2 

— oc m , 



(15) 



whereas for efficient ones it is proportional to rh. We can 
thus infer a relationship between the internal radius of the 
standard accretion disc (normalised to R g ) and the bolomet- 
ric luminosity by combining Eqs. 1141 and [TS1 



(16) 



in the inefficient case, whereas it goes as (,Ri n /R g )- 3/2 when 
the flow is efficiently radiative. 

The normalisation constant is actually undetermined in 
the relationship linking the radius and the typical frequency 
(Eq. I13p , and is subject to variations in the other equations 
used. Hence we plotted the dependencies demonstrated in 
equation [16] on Fig. 1161 using either the hypothesis that the 
frequency v\ used in Eq.[T2]is exactly equal to the Keplerian 
frequency (Eq. 1131 dashed line in Fig. 11611, or t o a sound 
wave. In this case, we have (see e.g. ICzerny||2004l . Eq. 3) : 



h d 

^sound — \ ] ^Keplerian- 



(17) 



This latter relationship implies that the normalisation 
in Eq. [16] is multiplied by a factor (hd/r) 2 . In Fig. 1161 
we plotted the corresponding relationship for (hd/r) =0.1 
(solid line), which seems to be more consistent with the 



points obtained by spectral fitting in the 10 
minosity range. 



10 -2 lu- 



6 DISCUSSIONS 
6.1 Summary 

In this paper we focused on the question of the evolution of 
the optically thick disc geometry during luminosity and state 
evolution observed in black hole binaries. Estimates of the 
internal radii values were obtained via modelling the thermal 
component in the spectra with a multicolour disc. In this 
simple model the luminosity is dominated by the internal 
part of the accretion disc that behaves like a blackbody of 
temperature Ti n and surface 4ivR 2 n . 

However, as the optically thick disc in a BHB is usually 
emitting in the soft X-rays or far UV, the normalisation of 
this component is strongly dependant the value of the in- 
terstellar (and internal) absorption along the line of sight. 
We studied four sources to determine if there were possible 
changes in the value of JVh during state transitions. Accord- 
ing to the data obtained with Swift, we demonstrated that 
when using a powerlaw for fitting the high energy compo- 
nent of the spectra, a decrease of the absorption is detected 
(above reasonable significance). This decreasing TVh is more- 
over observed when the source moves towards harder state. 



This effect is however not signi ficant when using the th ermal 
comptonisation model comptt jTitarchuk et al.lll994 ). 

Taking this trend into account, we then studied the evo- 
lution of the disc geometry, firstly in the particular case of 
XTE J1817-330. We demonstrated that this source was in- 
deed exhibiting t he beginning of a disc recession . This was 
already noted bv lGierliriski, Done fc Page! (|2008t ). but with- 
out taking into account disc irradiation, and contrary to 
what was stated by R07. This turnover in the disc behaviour 
is statistically significant (see Tab. [7}, even when fixing the 
value of the hydrogen column density. 

We then extended our study to a sample of six black 
hole binaries and tried to draw conclusions on the overall 
trend of the disc geometry evolution. These data, from a 
variety of instruments, are still consistent with a recession 
of the inner disc radius when the bolometric luminosity and 
the inner temperature decreases (see Fig. I10l and ll2[) or the 
spectrum hardens (Fig. I13p . 

Moreover in order to link the quiescent and high lumi- 
nosity data, we also demonstrated that Lboi had to scale 
with the square of the accretion rate flowing through the 
optically thick disc rh\ (Fig. 1 14p . Finally this evolution is 
consistent with accretion models where evaporation of the 
inner part of the disc is taken into account, however the 
strong ADAF hypothesis cannot be ruled out. The evolu- 
tion of the disc radii values in function of the bolometric 
luminosity is also consistent with the expected behaviour 
derived from timing in the lowest accretion states (see Eq. 
I16[) . We therefore propose in Fig. [17] a general sketch dis- 
playing the evolution of BHB in function of their state and 
accretion rate. 



6.2 Known limitations of the analysis 

We only focused our analysis on the strength of the disc 
component in order to determine the disc geometry. Future 
investigations should correlate these results with those ob- 
tained either by iron line width fitting (but see hereafter the 
known limitations) or studies on the timing characteristics. 

More importantly, we cannot exclude the possibility 
that the inner radii values obtained by using the simplistic 
multicolo ur disc model may so metimes be unphysical. Fol- 
lowing e.g lSobczak et all (|2000T l. a sudden flare could result 
in an apparent decrease in the inner disc radius value when 
using the simple multicolour disc to fit the thermal compo- 
nent when no such decrease takes place in reality. This could 
be due to increased spectral hardening and/ or Compton up- 
scattering of soft disk photons, thus depleting the thermal 
component, and hence decreasing the observed radius value. 

However, an increase in the radius when the source 
fades as observed in our study, seems more difficult to inter- 
pret in that framework. When moving towards the harder 
and lower states, one should on the contrary expect a de- 
crease of the radius value, according to the previous rea- 
sons. Therefore if we had taken those effects into account, 
we would have even obtained bigger radii value during the 
tra nsitions. This would lead t o equi valent results as obtained 
by iGierlihski. Done" fc Page! (|2008h . who took into account 
the irradiation of the inner part of the disc by hard X-rays 
coming from the corona on the same XTE J1817-330 Swift 
data than ours. They concluded that instead of receding by 
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Figure 17. Expected behaviour of a BHB during outburst in the DFLD. Note that this general shape will not be valid if the source 
transit by a very high state (or "Steep Power Law") with a high contribution of the powerlaw component. 



a factor two or three, as we obtained in this work, the inner 
radius would increase by a factor 6 to 8. 

This irradiation model could therefore influence the 
relationship between Lboi and R ln , and hence, the deter- 
mination of the accretion models performed in Sect. 15.2.11 
via a flat tening of the actual relation. This model de- 
scribed in iGierlihski. Done fc Page] (|2008T ) depends mainly 
on one extra parameter, the ratio R c = L c /Ldisc {L c be- 
ing the comptonised component luminosity), which deter- 
mi nes the deviation from the mu lticolour disc (see Fig. 6 
of IGierlihski. Done fc Page! |200Sj ). The higher the R c the 
more the deviation. As R c can be linked to the value of 
the powerlaw/disc fraction PLF (L c ~ -Lpl) giving R c ~ 
PLF /{I — PLF), we can perform a rough estimate of the 
influence of irradiation on Ri n in our study via the DFLD 
plotted in Fig. [T5] 

For XTE J1817-330, as in observations 15, 16 and 
17, PLF ~ 0.5 and therefore Rc = 0.5. Using Fig. 6 of 
IGierlihski. Done fc Page] (|200&f ). the radius should increase 
by a factor 1.25 only. This correction factor is negligible for 
points in higher states. The effect is however more striking 
for GX 339-4, SWIFT J1753.5-0127 and XTE J1118+480 in 
outburst. Indeed, both T08 data points have a R c ~ 4 and 
thus the radii obtained should then increased by a factor 3, 
giving R in ~ 66 and 42 R g . For SWIFT J1753.5-0.127, as 
R c can reach values of 5, it would imply a correction fac- 



tor of about 4. The case of XTE J1118+480 data obtained 
with Beppo-SAX is even more characteristic, as the lumi- 
nosity is completely dominated by the optically thin com- 
ponent and R c ~ 20. A s the l ack of coverage in Fig. 6 of 
IGierlihski. Done fc Page] (|200ct ) does not allow us to infer 
the correction in radius, this implies, an increase by a factor 
five should be a minimum. Paradoxically, this could solve 
the discrepancy observed between radius values determined 
from Beppo-SAX and EUVE data. The EUVE data indicate 
that the correction factor in radius should not be higher than 
2 (Rc ~ 2). On the contrary, both quiescent points might 
be slightly affected by irradiation as in that case R c < 0.2. 
In any case, our study could be improved in the future by 
investigating the effect o f adding irradia tion, or reflection 
models to the fits (see e.g lReis et al. I l2008l for a reanalysis of 
M06 data with reflection models, and iHiemstra et al. I l2008l 
for the case of SWIFT J1753.5-0127 in the lowest states). 

In addition to our extensive study, we note that an 
apparent disc recession during the decline was also noted 
in RXTE data at th e end of the observations analysed by 
ISobczak et"afl (|2000l ) on XTE J1550-564 (the inner radius 
value rises by a factor three when the total flux drops by a 
factor 50). 

The value of Ri n obtained in quiescence can also be 
subject to systematic deviation, but for other reasons than 
those mentioned so far. As already pointed out by (MC03a), 
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the optical/UV excess observed in quiescence can also be 
mimicked by an ADAF, without requiring a disc, if w e use 
the latest version of this model (jNaravan et al.iri997l ). The 
value of the transition radius R tI obtained by MC03a is then 
~ 10 4 i?s- As those latter authors do, we can hence reason- 
ably estimate that in quiescence 2 x 10 3 ^ 7?i n /R g ^ 2x 10 4 . 

Another problem arose from our study concerns the 
discrepancy observed in the outburst of XTE J1118+480 
between Beppo-SAX and EUVE data. We already men- 
tioned that disc irradiation could be one explanation, how- 
ever it does not seem to be sufficient to account for the 
very large differences observed. Another possible explana- 
tion could come from the uncertainties in evaluating the 
extinction in EUVE data and thus the EUVE slope (see 
e.g Fig.3 of (MCOf)), which could strongly affect the flux 
emitted by the source and hence the estimate of the radius. 
Independently to the discussions presented in (MC01) and 
(Ch03), we note that such a high radius from the EUVE data 
implies certain issues regarding the accretion efficiency. Fig. 
1 151 shows that in order to account for the observed data, the 
accretion rate flowing through the disc must be quite high 
(rha ~ 0.1). However, the low luminosity observed would 
also imply very inefficient flow, which is quite unusual com- 
pared to what is observed in the other BH binaries for such 
accretion rate. 

We note however that the case of XTE 11 18+480 in 
its 20 00 outburst has recently been revisited bv lReis et al.l 
(|2009h . These authors base their analysis on Chandra and 
RXTE data alone and conclude that a non recessed disc is 
present, attributing the soft excess observed in EUV and 
radio to the synchrotron component of the jet. However the 
handling of the absorption in this study seems puzzling. In 
the final fits, the value of the absorption seems to stick at 
its higher boundary (fixed at 1.3 x 10 20 cm -2 ). We note 
furthermore that when it is allowed to vary more freely, its 
final value is quite higher (~ 2.7 x 10 20 cm" 2 ). 



6.3 Consequences of a truncated/non truncated 
disc 

Assuming we state that the overall trend of the optically 
thick disc is to exhibit a progressive recession of its inner 
radius, how can we interpret the recen t observations of G X 
339-4 (M06) and SWIFT J1753. 5-0127 (|Miller et al.ll2006bh ? 
First, as stated by those authors themselves, it can be pecu- 
liar case, or perhaps observations performed when the source 
was not steady. We cannot exclude, due to the long dura- 
tion of the observation that the state could be a "mix" of 
a source transiting from the low to the high state for GX 
339-4 for example. 

Regarding the GX 339-4 observations from (M06) in 
detail, several points seems interesting to investigate. As we 
demonstrated (in Fig. 1 10 p , the bolometric luminosity is high 
(~ 0.1 I/ Edd) though with a low disc co ntribution. Hence, 
following iGierlihski. Done fc Pastel ()2008h . a complete treat- 
ment with the irradiation taken into account could give quite 
higher value to the inner radius observed as the powerlaw 
flux is completely dominating (R c ~ 112). Independently 
from this, the results obtained by the iron line fit with the 
laor model gives a very low inclination angle, i = 18°. Thus, 
if the inner disc and the orbit of the secondary are coplanar, 
as the mass function is ju = Afi sin 3 (i)/(l + M2/M1) 2 = 



5.8 ± 0.5 jHvnes et al.ll2003l . where Mi: BH mass and M 2 : 
secondary mass), this would imply a black hole mass of 
Mi ~ 2OOM (for M 2 = 1M ) or Mi ~ 215M Q (for 
M2 = IOMq). These values are quite unusual for an X- 
ray transient, and the only way to disentangle this is to 
sup pose the presence of a warped d isc in GX 339-4 (see 
e.g iFragile. Mathews fc Wilson 2001 for possible evidences 
in GRO J1655-40, or iMaccaronei 120021 for a more extensive 

study). 

Th e case of SWIFT J1753. 5-0127 in iMiller et al.l 
(|2006bl ) is more subtle, even if as in GX 339-4, the pow- 
erlaw flux is dominating (R c ~ 18) and thus the irradiation 
could have an effect on th e evaluation o f _Rj n . H owever, if the 
luminosity is quite low in IMiller et al l (|2006bl ). it mostly re- 
lies on the values of the mass and the distance of the source, 
which are currently not well known. If, we choose the values 
of our study (see Tab.[T}, we obtain Lboi/LEdd ~ 1-5 x 10 -3 . 
However, if d = 8.5 kpc, Lboi is nearly an order of magnitude 
higher. 

We can also note that the addi tion of a reflect ion 
component on the disc (such as in iReis et al ] 120081 or 
iHiemstra et al.ll2008l ) could also explain part of the soft X- 
ray emission, and thus could decrease the disc intensity and 
the inner radius value inferred from the fits. 

More fundamentally, a non truncated disc in the 
lowest state could raise som e issues according to 
iDubus. Hameurv fc Lasotal (|200ll ) simulations. When a disc 
is present near the ISCO, instabilities in the disc seems to 
be unavoidable and a truncated disc prevents the triggering 
of such outbursts. 

Finally, if we extend this study to larger scales, some 
data can also be obtained in the gap between 10 -3 and 
10 -4 I/Edd thanks to X-ray-Bright Optically Normal Galaxy 
(XBONG) and lo w luminosity AGN (M81 NGC 4579) data. 
Indeed, following lYuan fc Narava 11 (2004), the correspond- 
ing inner radius obtained are close to R ln ~ 10 2 R g at 
Lhoi ~ 10~ 4 I/Edd- Of importance, Fig. 3 of lYuan fc Naravanl 
|2004l ) seems indeed to exhibit roughly the same Lboi vs Ri n 
relationship we obtain in Fig. 1161 the corresponding "quies- 
cent AGN" being Sgr A* and M87. 



7 CONCLUSIONS 

Based on several datasets obtained with instruments sensi- 
tive at low X-ray-energies (SAX, Swift, EUVE), we studied 
the optically thick component in BHB spectra and its evolu- 
tion during outbursts via its contribution in the soft X-ray 
spectra. Any possible degeneracy with other processes such 
as the absorption in the line of sight was analysed. Depend- 
ing on the nature of the high energy emission process, a 
possible variation of the hydrogen column density during 
outbursts is detected. 

As a result of the good coverage by Swift, the case of XTE 
J1817-330 is examined in detail. Whereas for high bolomet- 
ric luminosities (> 10 -2 LEdd) the inferred inner radius re- 
mains fairly constant, a sudden increase of its value is signif- 
icantly detected when the luminosity is low (< 10~ 2 LEdd)- 
When including results coming from other sources and in 
quiescence as well, the hypothesis of a beginning of disc re- 



cession under 10 



10 LEdd is still favoured. The ob- 
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served disc recession rate with luminosity is also consistent 
with X-ray timing behaviour of BHs. 
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Table A. Table of the relevant fit/literature results we used in this study. The column densities values are expressed in 10 22 cm -2 , kT ln 
in keV, iodise in 10 3 X (Ri nj km/dio kpc) 2 cos(i), and the fluxes in terms of 10 — 10 erg. cm - 2 .s _1 . The powerlaw flux values in parenthesis 
were not taken into account as the photon indices values obtained are irrealistic. Absorbed powerlaw fluxes are computed between 0.05 
and 200 keV where unabsorbed disc fluxes are computed on the whole electromagnetic spectrum. Values written in italics are fixed during 
the fitting process. 
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00030919002 


0.40612 


1 18 1 
(117 





504+0. 018 
dJ -0.022 


3 


44+0-53 
^-0.44 


92.6+4I 


(0 


71+1.75-, 
"--1.86J 


(169-44l? 6 ° 9 ° 44 ) 


301 


,6/278 


Swift- 


XRT 


00030919002 


0.43712 


082 
027 





601 +0.018 
DU± -0.017 


3. 


1n +0.38 
lu -0.33 


87.7l? 2 & 8 


2 


9K+0.82 
ZD -1.70 


8.67l°;« 


305 


,8/285 


Swift- 


XRT 


00030919003 


0.45712' 


047 
(120 





rrofi + O.Oll 
oyD -0.011 


2 


1Q +0.16 
lu -0.17 


57.2±H 


2 


qq + 0.36 

,3 ' ;> -0.42 


90+0.00 

1*.ZJ_3 47 


439 


,5/425 


Swift- 


XRT 


00030919003 


0.46112 


033 
(125 





cm +0.009 
Ja -0.009 


2 


1S +0.14 
18 -0.15 


57.5l* : J 


2 


qq + 0.24 


17.70l2;f 6 


561 


,4/480 


Swift- 


XRT 


00030919004 


0.50412 


04(1 
043 





505+OOIO 
Ja,J -0.009 


1, 


66 +0.19 
DD -0.24 


45.01" 


2 


66 +0.21 
,DD -0.27 


17.46l° 2 ;|| 


600 


,5/494 


Swift- 


XRT 


00030919005 


0.43412 


150 
(144 





55Q+0.017 
dJa -0.017 


1 


cn + 0.57 
DZ -1.71 


97.7l2 7 9 6 


2 


5Q+0.78 
OU -0.52 


24 qq+5-49 
z<±.oa_ 24 39 


361 


,4/323 


Swift- 


XRT 


00030919006 


0.42112 


072 

035 





5 47 +o.()24 
d4 ' -0.023 


2 


77+0.43 
' '-0.38 


53.7+6-? 


2 


.9+0.59 
,lz -0.76 


19.69lf 9 2 3„ 


317 


,6/293 


Swift- 


XRT 


00030919006 


0.46112 


088 
057 





54 o+0.010 
dyi -0.018 


2 


gq+0.43 
-0.60 


50.4l 8 9 j 


2 


57+0.43 
°'-0.56 


16.98iJ;™ 


352 


,3/324 


Swift- 


XRT 


00030919007 


0.41112 


088 
032 





441 +0.032 
44 -0.034 


2 


7Q+0.98 
la -0.65 


22.913° 


1 


qo+0.28 
,au -0.31 


01 00 +8. 79 


281 


,8/291 


Swift- 


XRT 


00030919007 


0.45712 


042 
(130 





q 7Q +0.037 
°' :: '-0.043 


5 


q5+3.03 
oo -1.60 


24.0l 4 ° 


2 


q K + 0.21 
,OU -0.23 


z 7-34l 2 ;^ 


314 


,1/305 


Swift- 


XRT 


00030919008 


0.42712 


024 
022 





qt-i +0.024 
OJ -0.025 


3 


qo+1.58 
3Z -1.02 


12.9±\l 


2 


Oq + 0.12 

,uo -0.12 


21.45+ 2 f 7 


402 


,7/415 


Swift- 


XRT 


00030919008 


0.44112 


022 
021 





qqi +0.023 
• 3 '-' -0.024 


5 


44+2-31 
**-1.41 


i4.oi;;°3 


2 


16 +o.io 

1D -0.10 


20.431};™ 


574 


,9/496 


Swift- 


XRT 


00030919009 


0.40612 


021 
021 












2 


q fi +0.06 

,oo -0.06 


22.9911:17 


256 


,3/241 


Swift- 


XRT 


00030919010 


0.35312 


02!) 
02(1 





q.q + 0.049 

o4,;> -0.046 


1 


oc+1.93 
8b -0.84 


5 fi +0.7 
°-°-1.8 


1 


70+0. 08 
,,u -0.09 


49.65+11 


487 


,6/464 


Swift- 


XRT 


00030919011 


0.34412 


051 
(140 





626 +0.144 
DZD -0.100 





ns +0.08 
U8 -0.04 


9 5+O.2 


(0 


Q7+0.33-, 
u '-0.48^ 


(1811.88_ lgll gg ) 


158 


,7/165 


Swift- 


XRT 


00030919011 


0. -131 12 


074 
062 





4r -o+0.093 
4d3 -0.079 





OO+0.33 
zz -0.12 


9 .+0.3 
z ' 1 -0.8 


(0 


cc +0.24\ 
■ 66 -0.28) 


(236.73l 2 ^;f 3 ) 


168 


,6/162 
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Table C. Continuing table [5] 



Telescope 


Obsld/Ref 




fcT in 


^disc 






r 












GX 339-4 (continuing) 












Cam ft YRT 


uuuouy iyuiz 


n okq+0.024 
U.OOO_q Q23 


- 


- 




i 

1 


SA +0.06 
' 84 -0.06 


1 R 9q + 2.75 


9/19 S /9^^ 
Z^Z.O/ ZOO 


Ctmft YRT 
O W IJL - ^\ 1\ 1 


uuuouyiyuio 


n cinc;+ 042 

U.OUJ_Q Q 3g 








1 


cc+0.09 
,DO -0.09 


1 5 9K+3-95 
oq vo~l"41.85 

zj.iy_ 11 A2 


Q1 9 /1 on 
yi.z/ iuu 


Qimft YRT 


uuuouyiyuio 


u.zuy_ 135 










4R+0-34 
' 4,:) -0.31 


9Q Q /I 

zo. y/ iu 


Swift-XRT 


00030919013 


o n7 +0.050 

U.OU/ _Q Q 4 g 









1 


g 4 +0.12 
,u ^— 0.11 


15.09jf.JS 


69.3/77 


Swift-XRT 


00030943002 


229+ ' 028 
u.zzy_ 026 








1 


07+O.O6 

■^'-o.oe 


u.3iti-y 


192.8/189 


Swift-XRT 


00030943002 


n 97 o+0.030 
u.zia_ 02g 








1 


A9 +0.07 
'* z -0.06 


12.32+^06 


160.8/184 


Swift-XRT 


00030953003 


n 9^+0.039 
u.zoj_ 035 








1 


dA+0.08 

°*-o.o8 


16.09+"? 


116.7/132 


Swift-XRT 


00030953007 


S05+ 043 

U.OUO_q 039 








1 


oq+0.08 
"^-0.08 


34.66l»;« 


99.5/132 


Swift-XRT 


00030953012 










1 


A o+0.06 
' 48 -0.06 


35.93lS;« 


185.9/184 


Swift+XTE 


Spec. 1 (T08) 


0.850_ ; 070 


n 17S+0-007 
U.l/»_ 007 


7fi +39 
'"-26 


27.3 


1 


gq+0.01 
,Dy -0.01 


27.00 


532.0/211 


Swift+XTE 


Spec. 2 (T08) 


0.900±g : iig 


n , .7+0.011 

u - io '-0.011 


qcr + 30 
00 -18 


4.6 


1 


R4+0-02 
D4 -0.02 


9.20 


280.1/211 


ASCA G+S 


42010010 


0.^3 


914+ 060 
u ' zi4 -0.049 


, 10+4.99 
1 " LO — 0.84 


51+ 012 


1 


RQ + 0.02 

' DO -0.03 


in qq+0.03 
1U.UJ_ 06 


1030.70/954 


ASCA G+S 


42010000 


0.43 


171 +0.042 

u -i ' 1-0.033 


7 7K +31.84 
'•' D -5.88 


1 44+0-00 
i -^ 4 -1.44 


1 


fio+0.01 
o,3 -0.02 


17.511^7 


1201.70/1174 


ASCA G+S 


43001000 


043 


0.279±S;ggg 


9 co+0.48 
z - b8 -0.38 


4 ,+0.20 


1 


fi7 +0.01 
D '-0.01 


46.33l»; 8 87 


2506.09/1469 


XMM+XTE 


M06 


n o 7 +0.05 


0.38±g;gi 


fi4+ - 08 
u - D4 -0.08 


2.9 


1 


46+0-01 
' 4D -0.01 


156.26 


3899.2/2256 



SAX L+M+P 
SAX L+M+P 
EUVE+XTE 
Idcm+CXC+SAX 
CXC+HST 


21173001 
211730012 
(MC01) 
(Ch03) 

(MC03a) 


nl 0+0.002 

U.U10_o 002 
ml +0.001 
U.Ull_o 001 

0.013 
nll +0.002 

U.Ull_o 002 
0.012 


XTE J1118+480 
0.06210-012 78 o.65t 3 6 3 5 5 3 / 2 

0.0661°;°°* 992.35±»?i;«| 
0.024j°;°° 2 7.c5 
0.024j°;°° 2 6.6j 2 ; 2 e5 
l.l±g;ae-3 5.4e6 


9 49+1-36 
z - 4z -0.77 
9 n fi +0.40 
z - aD -0.86 

51.91 
44.85 
1.72e-3 


1 7/10+0-007 
l-'4o_ 007 
1 7R9+0-006 
1- ,DZ -0.007 
1 7S9+0-005 
1- '° z -0.005 

1 s+o-i 

2 02+0- 16 
z - uz -0.16 


07 1 O + 0.52 

°'"" J -o.69 
28.731";^ 

41.7 

40. 

1.8e-4 


213.65/189 
237.99/189 


CXC+HST 


(MC03a) 


0.019 


7.7e-4 


A0620-00 

2.8e8 


2.3e-2 


2 26+ 18 
z ' zu -0.18 


4.3e-4 





